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SUMMARY 


A study was performed to continue development of computational techniques initiated 
under Contract NAS9-10956, "Computational Techniques for Design Optimization of 
Thermal Protection Systems for the Space Shuttle Vehicle, " to be used in the studies 
for the Space Shuttle Thermal Protection System. The resulting computer code was 
then used to perform some additional optimization studies on several TPS configurations. 
The program was developed in Fortran IV for the CDC 6400, and it was converted to 
Fortran V to be used for the Univac 1108. Documentation of the study consists of the 
user's manual and the final report. 


The computational methodology is developed in modular fashion to facilitate changes 
and updating of the techniques and to allow overlaying the computer code to fit into 
approximately 131, 000 octal words of core storage. The program logic involves 
subroutines which handle input and output of information between computer and user, 
thermodynamic, stress, dynamic, and weight/estimate analyses of a variety of panel 
configurations. These include metallic, ablative, RSI (with and without an underlying 
phase change material), and a thermodynamic analysis only of carbon-carbon systems 
applied to the leading edge and flat cover panels. Two different thermodynamic analyses 
are used. The first is a two-dimensional, explicit procedure with variable time steps 
which is used to describe the behavior of metallic and carbon-carbon leading edges. 

The second is a one-dimensional implicit technique used to predict temperature in the 
charring ablator and the non-charring RSI. The latter analysis is performed simply 
by suppressing the chemical reactions and pyrolysis of the TPS material. 

Two types of stress analysis are also performed. A discrete element technique 
describes all of the configurations (with the exception of the carbon-carbon leading 
edge) and a more rigorous finite element analysis is used to solve the case for the 
large deflections of the RSI and its accompanying glue line. A dynamic routine 
predicts not only the noise excitation due to a number of external sources and the 
fatigue life of a number of simple panels but also a panel flutter analyses for the same 
simplified configuration. The weights of all configurations were determined by a 
unique new method in which total system weight per unit area is determined as the 
sum of the weights of component parts. Two cost analyses are provided: in the first, 
manufacturing costs are prescribed as the sum of the costs needed to fabricate and 
assemble all of the component parts utilized in the weight model. In the second 
technique, the total program costs are predicted based on historical data for simpler 
systems . 

Two basic types of input data are provided based on trajectory data: the first includes 
vehicle attitude from which aerodynamic heating rates and pressures are predicted. 

In the other heating rates and pressures are given as a function of time. The former 
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values are given in both hot wall and cold wall terms as options for the ablation 
routine. Standard program output includes heating rates, temperature, and stresses 
for discrete elements of the TPS analyzed as well as dynamic stresses and the panel 
weight and cost. 

Optimization and sensitivity studies are performed by the user by varying panel size, 
material properties, and configuration in a series of computer runs. The program 
sizes panel and insulation thicknesses. An optimum design is then identified as the 
one giving either minimum weight or cost as a function of the parameters being varied 
for the investigation. Sensitivity studies are performed by noting the change in system 
weight or cost due to the variation in some independent variable such as trajectory or 
heating prediction method for an optimum panel configuration. 
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SECTION 1 


INTRODUCTION 


The study performed under this contract was a follow-on to that accomplished on 
Contract NAS9-10956, "Computational Techniques for Design Optimization of Thermal 
Protection Systems for the Space Shuttle Vehicle." The output of the former study was 
a computer program which sizes a number of different configurations of metallic 
re- radiative thermal protection system (TPS) panels and supports to stipulated 
constraints of stress values, sub-structure allowable temperatures, and acoustic 
fatigue, among others. The development methodology is reported in Reference 1, and 
the detailed program description and user's instructions in Reference 2. 

As a result of the continuing space shuttle systems studies pursued by NASA and 
Aerospace Industry Contractors, the interest in metallic re-radiative TPS waned, 
and was replaced by concepts employing carbon-carbon, Reusable Surface Insulation 
(RSI), or charring ablators. The follow-on study concentrated on these latter concepts, 
and also included study and implementation of improvements to the basic computer 
program to reduce computational time and to incorporate the latest aerothermodynamic 
methods recommended by the Thermo Panel of the Space Shuttle Aerothermodynamic/ 
Configurations Working Group, noted in Reference 16. 

Portions of the development of the analytic methods, equations and data banks were 
generated under Company sponsored activities. These were made available to be 
programmed and incorporated into the TPS computer program. The results are 
included in this report because of their applicability and interest to NASA. 

The study was performed in a series of related tasks. These investigations included: 

(1) the improvement of existing analyses already developed for metallic thermal 
protection system cover panels and the development of four subroutines for treating 
passive and active TPS's, including carbon-carbon composites, RSI, charring ablators, 
and a heat sink substrate using phase change material (PCM), (2) computer program 
applications which include both short term improvements to speed the computer use 
and improve its efficiency as well as investigations to assess the desirability of 
incorporating formalized optimization procedures into the system sizing routines, 

(3) performance of preliminary optimization and sensitivity studies to demonstrate 
use of the techniques, and (4) the documentation of the results of the studies and the 
improved computer program. 
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SECTION 2 


COMPUTER PROGRAM REFINEMENTS 


The prediction methods and analyses developed for use in the TPS sizing computer 
program as part of this year's effort are described and discussed in this section. 

The complete computer program, including those sections developed under last 
year's effort which remain unchanged, is presented in a companion volume entitled 
"User's Manual" which contains the listing and operating instructions. 

2.1 MATHEMATICAL AND NUMERICAL METHODS 

2.1.1 THERMODYNAMIC ANALYSIS . For the case of the metallic cover panel, heat 
transfer is determined two-dimensionally. The temperature at any point in the 
structure is a function of external and internal convection and radiation rates and 

the conductivity and thermal inertia properties of the structure itself. In this program, 
structural temperature distributions are evaluated through use of the lumped parameter 
method of finite differences. 


The surface and/or structure is divided into an arbitrary number of small segments. 
The segments are arranged in rows parallel and columns perpendicular to the surface 
as shown in Figure 2-1. 


AERODYNAMICALLY HEATED 
OUTER SURFACE 



Figure 2-1. Surface/Structure Segmentation 


Then, for some small time increment (At), the net heat flux to each surface segment 
is determined. For a segment i, then, the temperature change from time to time, 
t + At, is 


T i (t+At)-T i (t) = 


q net A BL" ^INSD" 2 ^ X *N " S(q rad A >N 


At 

(WC p ). 


( 2 - 1 ) 
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where 


A^j = area perpendicular to direction N 

k N = effective thermal conductivity in direction N 

Xjq = length of conduction path in direction N 

ATn = temperature difference between adjacent elements in direction N at time t 
cirad = heat transfer by radiation between nodes 
W = weight of element i 
Cp = specific heat of element i 

It is assumed that all the mass is concentrated in a point at the centroid of the segment, 
An is equal to the segment interface area, and Xn is equal to the distance between 
centroids in direction N. 

The net heat transfer is given by 

q net %ons ^i - q rad (2-2) 

where q* is the boundary layer convective heat transfer rate, q^d represents the energy 
loss due to surface radiation, and q cons is an optional multiplying factor. This factor 
may be used to approximately allow for the effects of shock wave interactions, flow 
divergence, etc. 

The term (qA)iNSD provides an optional capacity to include internal convective cooling 
of the structure. The cooling occurs on the backface of the last segment. 

Terms in Equation 2-1 which do not apply to a given element are dropped out for that 
element. Thus, the q ne t applies only to elements of the first row, i.e. .those representing 
the surface, and the (qA)jNSD term applies only to elements of the last row, i.e., those 
representing the backface. Internal radiation heat transfer is also taken into account 
in this program. The accuracy of Equation 2-1 is dependent upon the size of the 
segments and the computation interval and improves as these parameters are decreased. 
The program may be used for either one-dimensional or two-dimensional arrays of 
segments. 

In the following, the two-dimensional case is discussed; the one-dimensional problem is 
treated as a special case of a single column. The structure is set up in a matrix shown 
in Figure 2-2. The material may not be homogeneous in either the x-direction or the 
y-direction. The maximum number of rows and columns is nine each. Configurations 
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(b) TYPICAL NODES 

Figure 2-2. Structure Heat Transfer Matrix 


are given in Figure 2-3, and nodal breakdowns for the thermodynamic and stress 
analyses are given in Figure 2-4. 

The temperature at any segment at time t + At is dependent upon the summation of 
direction heat transfer rates at time t. The general equation is given by Equation 2-1. 

From Figure 2-2, for a typical node, the temperature, T, changes from time t to t+ At 
by 

T 1J (t+At)«T lJ (t)+^ -gjOK (2-3) 

Ajj is the thermal mass of the node as defined in Table 2-1, together with other parameters. 
The equations for computing Qk's are listed in the following: 

a. Conduction heat transfer 

= < T i-l, j - T i, j> / < 2 ‘ 4 > 
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2. FLAT CORRUGATION WITH SKIN 



— -J h-TC 


3. RIB-STIFFENED PANEL 



4. SKIN-STRINGER 



5° 0PEN CORRUGATION 6. OPEN CORRUGATION 

(CIRCULAR ARC CORRUGATION) 


Figure 2-3. Panel Geometries 
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Figure 2-4 . Configurations for Thermodynamic and Stress Analysis 
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Table 2-1. Notations for Heat Transfer Calculations 


Symbol 

Calculated Term 

x! 

x i 

X 

yl . 
l.J 

y j 

A i,J 

(p Vi.j x j y j 

B i»j 

x i / < 2K i,j y j> 


Symbol 

Calculated Term 

c i,i 

V (2K i.fi* 

D u 

- 

e i,j 

B i,j +B i+l,j 

f m 

C i» j + C i, j+1 


Qs'fTM-l-Tyt/Pi.j.! 

Q4'(Ti,j + l-T 1J )/F lij 

b . Radiation heat transfer 


Q K = a EL 3 . ( T 4 

a > D m,n 



where 


L = x| (for k = 3 , 4) 

= yjj (for k = 1,2) 

c. Nodes on the first column 


(2-5) 

(2-6) 

(2-7) 


(2-8) 


Q 1 = 0 

d. Nodes on the last column 

Q 2 = 0 

e. Nodes on the first row 

q 3 = *1 %eti 

f. Nodes on the last row 

Q 4 = h ins x i ( T ins " T ij) 
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One particular improvement which was made to the heat transfer calculations of the 
metallic TPS's was the adoption of a time step variation. Since the numerical 
technique for solution is explicit, small time steps must be taken at various times 
during the trajectory to avoid numerical instability. However, at other times, 
conditions are changing so slowly that small time step sizes are prohibitively 
expensive. The temperatures of each node are now monitored throughout the calculations, 
and when it is determined the largest temperature change is small, the time step size is 
lengthened still however being maintained within a theoretically established stability 
criterion. The converse is also true. When temperatures change rapidly, the step 
size is decreased. Such variations lead to a 39% savings in computer time with no 
noticeable loss of temperature accuracy. 

For the case of the thermodynamic analysis of the ablator and RSI TPS's, the one- 
dimensional heat transfer computations are performed by an implicit technique which 
avoids the difficulties of numerical instability. Since this technique was adopted almost 
in total from an existing computer program, the reader is referred to the original 
reference by Curry for details (Ref. 3). 

2.1.2 STRUCTURAL ANALYSIS . During early stages of the contract, a separate 
structural analysis was performed on the supporting heat posts which attach the simply 
supported cover panels to the primary structure of the space shuttle vehicle. Results 
of a parametric study of various loading conditions described below, demonstrated that 
for the flexural rigidity established by acoustic test, the design point for this part of 
the structure is minimum gauge. Hence, for the remainder of the effort, these values 
were used. 

Analyses of Supporting Structure 

Sizing of the tubular supports for both static and dynamic loads of the existing metallic 
TPS support configurations is not performed by Program P5490. However, an existing 
computer subroutine was adapted to conduct this analysis. As is the case for the TPS 
panel, the posts are sized for both static strength for aerodynamic loads and dynamic 
response due to the acoustic environment. The vibrations induced by the engines, which 
are transmitted mechanically through the shuttle structure, are omitted in the analysis 
because the excitations are attenuated rapidly except in the immediate area of the thrust 
structure. The dynamic response of the post, however, to the acoustic environment is 
important because it has been demonstrated at Convair Aerospace in the design of TPS 
supports that the dynamic stress can be a critical item (Reference 4). 

The design requirements applied to this analysis are: 

1. The applied ultimate 3-a stress on the tube wall must not exceed a specified cut- 
off value (examples are 78, 000 psi for TDNiCr and 120,000 psi for titanium). 

2. The applied 3-a stress must not exceed the crippling stress of the tube walls. 
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3. The applied limit 1-c stress must not exceed the allowable fatigue stress 
determined from the random S-N fatigue curve for an input required panel 
life. 

4. The post must have a critical column buckling load greater than the applied 
ultimate axial load acting alone. 

5. The post wall thickness must be greater than input minimum gauges. 

The input data to this particular post sizing analysis includes 

post length 

weight of the dynamic forcing function 
pressure force 
ultimate factor 
S - the sigma factor 
K - buckling constant 

vibration environment (power spectral density as a function of frequency) 

S-N data for the post material 

The heat shield support posts are represented as a single degree of freedom system 
whose natural frequency is idealized as that of a massless uniform cantilever beam 
with a concentrated mass at the free end (Figure 2-5). The loads on the post include 
an applied axial load P and a lateral load Q rms . The former is due to the differential 
pressure between the external aerodynamic load and the internal pressure behind the 
TPS. The applied lateral load W (transferred to the TPS panel and post by the acoustic 
environment) is a function of panel weight subject to a specified rms acceleration at 
the system's fundamental frequency 

f n = ~ ^ (1158 EI)/(L3W) (2-9) 

A power spectral density curve specified to describe the environment (a typical one is 
illustrated in Figure 2-6) then gives the rms acceleration as 

g rm B -/f fa f (2 _ 10) 

where 

E = modulus of elasticity 

I = moment of inertia 

L = tube length 
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Figure 2-5. Support Post Structures 
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Figure 2-6. Typical Interface Liftoff Vibration Specification 
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S n = amplitude of the power spectral density function in g 2 /cps at the 
frequency f n 

Q = damping ratio 

The rms lateral load then becomes simply 

Qrms = ± Srms ‘ W ( 2-1 l) 

In order to size the uniform tubular post at its base, the rms bending moment is given 
by beam column analysis as 


M = QJ tan (L/J) 


( 2 - 12 ) 


where 


Q 

J 

V 


n 


= lateral load 
= j/( 7}EI)/(UF p") 

= plasticity correction factor 


/ 


i + l«<w s o >“' 1 


= slope of the plastic stress-strain curve 


Smax = specified cut-off value for maximum stress 
S Q = specified base stress 


(2-13) 


UF = ultimate factor (e.g. , 1.4 to 1.5) 


The maximum applied stress a mny for the applied 3-a stress is given by 


P + S 


4MD 


irD®t 


max Dt 

where S = 3 for the 3-a applied stress 


(2-14) 


Thus Equations 2-12 and 2-14 are implicit in the plasticity correction factor. The 
performance of this calculation necessitates iteration. The 1-a alternating stress 
to be used in the fatigue analysis is given simply by 


a 


alt 


8M D/2 
itD 3 t 


(2-15) 
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The allowable crippling stress is given by 


F cr = K 7j E (t/D) (2-16) 

where K = buckling constant, and the critical column load is given by the Euler equation, 
modified for plasticity 


„ _ C 77 7) El 

cr 

L 2 


(2-17) 


where C = fixity coefficient, which for a column fixed at one end and free at the other is 
0.25. 


The design requirement for the latter load is that the critical column load exceed the 
product of the applied axial load and the ultimate factor. For this particular analysis, 
the allowable random stress Of a t corresponding to the 1-a applied stress is computed 
from an allowable S-N curve (Figure 2-7 is an example) for the number of stress 
reversals given by 


N = 4 N f T f n (2-18) 

where 

Np = number of flights 

T = duration of the vibration environment per flight and the factor 4 is a 
factor of safety for fatigue life scatter. 

The design factors corresponding to design requirements 1 through 5 are then computed 
as 

Design Requirement Design Factor 

* CT max^®max 

2 a max/ F cr 

3 a alt/ CT fat 

4 (UF P)/P cr 

5 t^/t 
where t^g is the minimum gauge thickness, an input parameter. 

The post design procedure of the existing computer program consists of iterating upon 
the post wall thickness t for a given tube diameter until the largest design factor is 
unity. All others must be less than unity, indicating that the applied load is less than 
the allowable , 1. 
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Figure 2-7. R = -1 and Random S-N Curves for TDNiCr 





Results of studies conducted for TDNiCr posts 
W =0.80 lb (vibrating mass) 

£ = 0.125 (damping ratio) 

P = 48 lb (compressive axial load) 

indicate that smaller post diameters yield lighter posts. This is shown by Figure 2-8. 

In addition, recent data developed at Convair have shown that this damping ratio may 
be low by an order of magnitude due to the frictional damping of this TPS of individual 
panels rubbing against each other. With this in mind, the entire system of cases 
represented in Figure 2-8 and a sensitivity study accompanying that investigation were 
re-run with a damping ratio of 10%. The results show that the design requirement in 
all cases run was minimum gauge; thus, the design philosophy of post supports alreay 
utilized in Contract NAS9-10956, to use minimum gauge posts, is confirmed. 

2. 1. 3 DYNAMICS ANALYSIS, PANEL FLUTTER AND STABILITY . Panel flutter is a 
self-excited instability of elastic panels in supersonic flow, during which the panel 
oscillates in a direction normal to its plane at a constant amplitude; flutter usually 
results in fatigue damage to the panel. Considerable theoretical and experimental 
research has been conducted to understand the phenomenon and to develop design criteria 
to prevent it (Reference 5). The design approach recommended by Lemley is incorporated 
into the TPS sizing routine. The minimum panel lliickness required to preclude panel 
flutter within the flight envelope is thereby established as a design constraint. Table 
2-2 lists significant parameters that affect panel flutter. 



TD NiCr POST DIAMETER (in.) 
Figure 2-8. Post Section and Weight vs Post Diameter 
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Table 2-2. Parameters Affecting Flutter Speed The flutter analysis is performed 

using many computations already 
made in the sonic fatigue analysis . 
The stability limit for a flat panel 
is given in Figure 2-9 with the 
panel flutter parameter $ plotted 
as a function of the panel length 
to width ratio. The curve is 
developed from experimental 
data (References 6 and 7). The 
function f(M) is given by Figure 
2-10. The design occurs at the 
maximum value of q/f (M) for the 
given trajectory. Once this has 
been established, the computa- 
tional procedure for establishing 
critical panel thickness is as 
follows: 

1. Determine $ for the panel length- to- width ratio from Figure 2-9. 

2. Determine q from flow field calculations. 


Flight 

Conditions 

Physical Data and 
Geometry 

Environmental 

Conditions 

Mach no. 

Young's modulus 

Inplane Stress 

Dynamic 

Pressure 

Length 

Differential 

Pressure 

Angle-of- 

Attack 

Width 

Differential 

Temperature 

Flow 

Angularity 

Length- to-Width Ratio 

Curvature 

Cavity 

Edge Conditions 
Thickness (to be 
determined) 



3. Determine f(M) from Figure 2-10. 

4. Evaluate the panel modulus of elasticity E from the panel material property 
values at the design point temperature. 


5. 


6 . 


Compute the panel thickness tg for flutter stability, where 


$ = 


" f(M)E ~l 1/3 

q 


|\B 


(2-19) 


Compare this value tg with the equivalent panel thickness already established 
in the fatigue analysis . 


7. If t B is less than the panel equivalent thickness, panel flutter presents no 
problem; if greater, the results are so noted in the fatigue analysis output, 
(Figure 2-11). 


2.1.4 WEIGHTS/COST ANALYSIS 

2 . l. 4. 1 Parts Listing. The parts listing procedure requires that a library be 
maintained to call out an associated list of detailed parts when a given panel or 
structure configuration is specified. Each part has associated with it a list of 
manufacturing processes and amounts of material required to produce the part. 
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MACH NUMBER CORRECTION FACTOR, f(M) 





















FPEOtJENCY 


577.9 H T 


POUNOARY LAYE° NOISE = 131.4 TP - REYNOLDS NO. = 3.4531E+P8, MACH NO. = 1.5315E+00 

OYNAMIC PRESSURE = 5. 0738 r -«- C 3 , VELOCITY = 1.4433E+03 

POCKE* ENGINE NOISE = 135. ft DP - APPA°ENT NOISE SOUPCE AT tl.O ET 

ARES NOISE = T7.6 OP 

JET SCPUPRING NOISE = 131,4 OP 

MAXIMUM PENDING STPfSS = 554 PSI 

MAXIMUM PMS O c FLECTION = . 00CC1 IN 

DOUNPARY LAYER NOISE CRITICAL STRfSS = 3.1053E+01 PSI 

NO. OF STRESS REVERSALS = 7.3451F+CE 

♦•♦PANEL IS GOOD FOF THIS CONDITION 

CRITICAL SIRES' MAY PE INCREASED TO 3.3714E+G4 PSI 

MAXIMUM PENDING S T RfSS = 933 PSI 

MAXIMUM RMS DEFLECTION = .OOCtl IN 

ROCKET ENGINE NOISE CRITICAL STRESS = 7.567AE+01 PST 

NO. OF STRESS REVERSALS = 1.9SC9F+CE 
♦♦♦PANEL IS GOOD FOF THIS CONDITION 

CRITICAL STRESS MAY pr INCPEAS.En TO 4.3595E+04 PSI 

MAXIMUM PENDING STRESS = 4 PST 

MAXIMUM PMS DEFLECTION = .00CC0 IN 

JET FL YPACK ENGINE NOISE cpttical stress = 3.0903E-C3 PSI 

NO. OF STRESS REVERSALS = 1.7997E4.07 

♦♦♦PANEL IS GOOD FOF THIS CONDITION 

CRITICAL STRESS MAY pr INCREASED TO 3.3849E+04 PSI 

MAXIMUM RENDING STRFSS = 1757 PSI 

MAXIMUM PMS DEFLECTION s .00003 IN 

JET SCPUPPINC- NOISE CRITICAL STRFSS = 3.3 859E + 03 PSI 

NO. OF STRESS REVERSALS = 7.1317F+C4 

*»»p ftN c L J5 GOOD fop THIS CONDITION 

CRITICAL STPESS MAY RE INCREASED TO 5.96R9E+C4 PSI 

COMPOSITE CRITICAL STPESS = 3.437AE+73 °SI 

NO. OF STRESS REVERSALS = 1.3393E+05 

♦•♦PANEL IS GOOD FOF THIS CONDITION 

CRITICAL STRFSS MAY RE TNCPEASFD TO 4.317PE+04 PSI 

PANEL fluTTEF ANALYST' - D'SIGN POINT AT 7.43Q0E+01 SECONDS OF FLIGHT 
YOUNGS MODULUS = ?.1600E+ n 7, MACH NO. = 1.1933E+0C 

0 LOCAL/FCM(AMl) = 1 • 0045E+ C 3 « Fhi = 4.0750E-G1 

♦♦♦PAMEL IS ST A «L r FO D FLUTTER 

PANEL EFFECTIVE THICKNESS = 3.3534E-01 IN. 

FANEL EFFECTIVE THICKNESS NECESSARY FOR FLUTTER STARILITY = 1.4853E-03 IN. 


Figure 2-11. Output of Fatigue and Panel Flutter Analyses 
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In this way, a method of deriving the manufacturing costs and material costs of the 
complete TPS system has been developed with its basis at the detailed parts level. 

The parts listing process occurs as a function of the panel concept (specified by 
KINDP-1, 2, 3 as shown in Figure 2-12). A detailed parts listing for the specified 
configuration is included in the manufacturing cost summary, Figure 2-13. 

2. 1.4.2 Weights . Actual and purchase weights are computed for each detail specified 
in the parts listing process; the results are listed in the output under the headings 
ACTUAL WEIGHT and MATERIAL WEIGHT. Actual weight is just what the name 
implies. It is computed based on the actual geometry of the finished detail part, taking 
into consideration all the necessary manufacturing and design requirements that normally 



RIB-STIFFENED HONEYCOMB SANDWICH CORRUGATED 


Figure 2-12. Panel Concepts 


THF9MAL COCT^CTTOK SY®TE»% ?PA G E SHUTTLE ST A XXXX 


rOKFISl'OATT^N PANEL T*PF 1, TORSUGATEC NOMINAL PANEL SIZE 3.0 X 2.0 FT 
STRUCTURE Type a NOMINAL r T A NT n FF LENGTH .5 IN 


THFro^TL ACTUAL MA T L S^O LA n 0R LABOP CV-HD M ATL « LAPGp OVERHO FACTORY MATERIAL F A PR IC A T 

peef-oxe CT WEIGH* WEIGHT WEIGHT HOUR? HOUR? PATE PATE PFR L° COST COST COST COST CCST 


FlNcL • 1 0.0110 tli « **7? 1 ", 71 9 0.0000 P.OO C.00 0.00 0.9P 0.0 0 0.0 0 0.0 0 0.0 0 0.0C 

CPPPUGN® 1 . o.oor 14.27? 17.914 9.6906 24.73 4 . 75 8.71 14.50 117.45 205.54 723.00 371.38 594,77 


CTPur*U B E 1 0.0P0 j , 14 P 1.701 o.nooc 0.00 P.OO 0.00 n.00 0.00 0.00 0.00 P.no 0.00 

OFAM IN A 1 O.P0P ,14® ,4?5 .338.1 .85 4.75 6.31 14.00 4.01 7.03 11.04 5.55 17,59 

op AM ru A 1 0,000 ,217 ,?") .2103 .53 4.75 8.31 14. 0C 2.50 4.77 6.87 4.32 11.19 

CORNER® 1 o.roo .05* .998 .0628 .16 4,75 8.31 l’.5C .75 1.30 2.05 1.30 3.35 

PC«T rCRN 2 0.000 , OF 7 .793 . 3111 .59 4,75 8.71 12.00 ?.79 4.87 7.65 5.05 12.70 

55® AL 1 0.PD0 .07? .11? . 05 P5 .13 4.75 8.3t 14.5P .50 1.05 1.65 1.79 3.44 

F A C T ENFRF 40 0.008 .775 .412 O.OOOC 0.00 T. 00 0.90 . 1’2.71 0.00 9.00 0. 00 60.22 60.22 


TOTAL THEORETICAL W®IGht 
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average material 
average l a°cr rate 
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224.16 * 

352. 2 B < 

3EC.61 * 

TP ?, 86 * 

PI. 05 « 

723.91 * SUB- A ?®FM9L Y COST 

U0.65 «rsn FT 


? . 57 

L9/SC FT 


STD 

TOTAL 

LA fi CP 

OV-HR 

LAPCR 

OVERHO 

ASSEMBLY 

18,73 

«/LE 

TA®N 

HOURS 

HOURS 

rate 

RATE 

COST 

GOST 

COST 

4.75 

8.31 

* /MO 
t/wp 

LOCATE 

.0456 

.11 

5.20 

1.75 

.59 

1.04 

1.63 

13.06 

* /HP 

CLAHF 
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.70 

5.20 

1.75 

3.64 

6.37 

10.01 

26,06 

«7MP 

PRILL 

.0416 

.IP 

5.20 

1.75 

.54 

.95 

1.49 

45.60 

8/LP 

SECURE 

.0570 

.14 

5.20 

1.75 

.74 

1.3U 

2.04 

.40 


INSFFCT 

.0304 

• 06 

5.20 

1.75 

.40 

.69 

1.09 

.40 


CISASSY 

, G R60 

.24 

5.20 

1.75 

1.25 

2.18 

3.43 

C.00 


CLEAN 

.0780 

.OR 

5.20 

1 .75 

.49 

.86 

1.36 


Figure 2-13. Manufacturing Cost Summary 
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go into producing a real part. Purchase weight is the weight of the raw material that 
must be purchased to produce each detail part. It is always larger than the actual 
weight. Calculation of the actual and purchase weights for a TPS detail part of the 
thermal protection system uses an equation of the form demonstrated below. 

Actual weight = density * length * width * thickness ( 2 - 20 ) 

Purchase weight = density * (length + Ci) * (width + C 2 ) 

* (thickness + C 3 ) (2-21) 

where C]_, C 2 , and C 3 are incremental distances added to the part geometry to account 
for the material removed during the cutting to size manufacturing operations. The 
actual equations can be found in the subroutine WTTPS in the program listing. 

2. 1.4. 3 Manufacturing Processes. To predict the manufacturing processes required 
for each detail part, a library of shop orders and shop planning records was established. 
These processes, along with the Convair Aerospace San Diego operation standard hour 
data, were studied and used to identify basic standard shop operations and to correlate 
each detail part in the parts list with a set of these standard shop operations. The data 
are used to compute the standard hours necessary for each shop process required in 
the production sequence of each detail part. 

Derivation of Standard Hour Equations 

The Standard Hour equations are derived based on Standards data acquired through the 
industrial engineering department. Typically the Standards data are comprised of three 
major items - a machine setup time, a part handling time, and a machine run time. 

These may be considered separately, or are sometimes combined into tables. Each 
item considered may be either a constant time increment per job or per part, or may 
be a function of the part geometry or characteristic machine operation parameters. 

For example, in edge burring, the handling time is a function of the part area in square 
inches, and the run time is a function of the length to be burred in inches. In straddle 
milling, the handling time is again a function of the part area, but run time is a function 
of the volume of material to be removed from each part in cubic inches. In each case 
the setup time is a constant used once per job. 

An example of the initial form of the Standards data is included in Figure 2-14. The data 
presented is in table form and represents the Standards for a HUFFORD A -12 extrusion 
stretch forming press, Convair machine code 8030. In this case the total standards 
hours are made up of two basic items, machine setup time (one increment per part for 

performing, and one for finish forming). The run time increments are a function of the 
overall part length. 

The development of the Standard Hour Subroutine involved acquiring the general Standards 
data and deriving an equation for each manufacturing operation based on the characteristics 
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SETUP: 0.52 (once per die change) MACHINE CODE: 8030 

PREFORM 

Length 
Std. Hr. 


1 

31-45 

.0315 

46-60 

.0345 

61-75 

.0375 

" 

91-105 
i 0435 

Length 
Std. Hr. 

106-120 

.0465 


136-150 

.0525 

151-165 
. 0555 

166-180 

.0585 

181-195 

.0615 

196-210 

.0645 

FINISH FORM 

Length 
Std. Hr. 

— 


31-45 

.0655 

46-60 

.0685 

61-75 

.0715 

1 

1 

91-105 

.0775 

Length 
Std. Hr. 

106-120 

.0805 

121-135 

.0835 

136-150 

.0865 

151-165 

.0895 

166-180 

.0925 

181-195 

.0955 

196-210 

.0985 

NOTE: Length in inches is based upon the bill of material length of part. All 
values include stock allowance for vice jaws. 


Figure 2-14. Standard Hours for Press, Extrusion Stretch 
Form 


process and part parameters. For the example Standards data, a general equation 
takes the form 

STDHR = 0.52 + N [ f x (L) + f 2 (L) ] (hours) (2-22) 

where 

0.52 = setup time (constant per job) 

N = total number of parts to be produced 

f 1 (L) = preform time as a function of part length 

f 2 (L) = finish form time as a function of part length 

The functions of length f^L) and f 2 (L) are determined by curve fitting the data in the 
Standards table. In this case a linear curve fit is sufficient and the functions 
resulting are . 

f 2 (L) = 0.002 L + 0.058 

f 2 (L) = 0.002 L + 0.024 

The resulting Standard Hour equation for this particular press forming operation is 
then 
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STDHR = 0.52 + N (0.004 L + 0.082) 


The Standards data are usually derived for aluminum only. To apply the data to 
additional materials, material complexity factors are utilized. The material 
complexity factors account for the difference in manufacturing time requirements for 
performing identical tasks or operations on different materials. These factors are 
typically required only for those manufacturing operations associated with material 
removal such as drilling, routing, burring, cutting, etc. Operations which usually 
do not require complexity factors are cleaning, layout, identification, painting, etc. 
Factors of 3. 8 and 4.2 are commonly used as complexity factors for steel and 
titanium respectively. Methods are available for deriving complexity factors for 
other materials; these methods will be explored and factors developed for all the 
program material selections. Figure 2-15 presents a list of the manufacturing 
operations currently available within the Standard Hour subroutine. Provision has 
been built into the subroutine to allow the future addition of any number of new 
operations. 


305 

Burring, Edge 

339 

Layout Holes, Sheet Metal 

308 

Clamping 

340 

Layout Holes, Machine Shop 

310 

Cleaning, Basket (Degrease) 

341 

Layout Part, Machine Shop 

311 

Cleanup, Assembly 

343 

Milling, Chemical 

312 

Disassembly 

348 

Milling, Straddle 

314 

Drilling, Single Spindle 

353 

Painting, One Coat 

315 

Drilling and Tapping or Reaming 

371 

Routing, Edge (Single Piece) 

316 

Drilling, Assembly 

375 

Sawing, Cross Cut (Extrustions, 

317 

Forming, Hot Pressure Forming 


Sections) 


(Sheet Metal) 

380 

Securing 

326 

Heat Treatment 

382 

Setup, Assembly 

332 

Identify, Rubber or Steel Gang Stamp 

383 

Shearing to Size 

336 

Inspection, Assembly 

385 

Surface Treatment 

337 

Inspection 

392 

Turning - Brazing 

338 

Layout Part, Sheet Metal 

397 

Welding - Brazing 


Figure 2-15. Basic Manufacturing Operations for Which S tandar d 
Hour Estimates are Available 

Seven standard shop operations were identified for this study 


Mfg. Operation 

Mfg. 

Mfg. Operation Mfg. 

Index 

Operation 

Index 

Operation 

1 

Cuttoff 

4 

Drilling, Routing, Deburring 

2 

Milling 

5 

Surface Treatments 

3 

Forming 

6 

Heat Treatments 



7 

Painting, Identification 
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Cutoff includes issue of the stock material to the shop and all cutting or sawing operations 
required during production. Milling includes milling, boring, and turning operations; 
forming includes all forming, stretch forming, and bending operations. Drilling, routing, 
and deburr ing operations are considered as a related group. Surface treatments include 
chemical milling, etching, anodizing, and peening. Heat treatments, painting, and part 
identification are also considered. An example of a callout for the manufacturing 
operations required to produce a corner post is 

Corner Post (Titanium) Manufac turing Operations 

KR and KS = 1 . 2 . 4 . 6 . 7 

Thus, the operations include cutoff, milling, drilling, heat treating and identification. 

The set of operation indexes KR and KS are then used by the program along with a 
material form index KK to call out values from the arrays KSETUP (KK, KS) and KRUN 
(KK, KR) located in the K-TABLES with BLOCK DATA. These values (in hours per 
pound actual weight) are then used to estimate a setup time and a run time required to 
manufacture each detail part. It is assumed that one basic machine setup is required 
per detail part, and that the setup for each additional like part is included with the run 
time. The equations take the following form 

Setup Time (SETUP): ACWT * KSETUP (KK, KS) 

RunTime (RUNTM): ACWT * KRUNTIM (KK, KR) 

where ACWT is the actual weight of the detail part. The setup and run times calculated 
are in standard hours as discussed in the following section. 

To predict the operations required to assemble the detail parts into the basic sub- 
components and then to develop a complete final assembly, a library of shop assembly 
planning records was established. Basic processes were identified and correlated with 
each subassembly and with the final assembly. This information is established in 
equations to compute the standard hours for subassembly and final assembly. 

Seven standard assembly operations were identified for this study. They are given in 
the tab le below along with their appropriate assembly operation. 


Assembly 

Assembly 

Assembly 

Assembly 

Operation Index 

Operation 

Operation Index 

Operation 

1 

Setup 

4 

Securing 

2 

Clamping 

5 

Inspection 

3 

Drilling 

6 

Disassembly 



7 

Cleaning 
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Setup includes the mounting of parts to assembly fixtures, and clamping involves the 
clamping together of the parts to be assembled. Drilling considers the location, 
center punching, and the drilling of all required holes. Securing consists of the 
actual attachment of required fasteners, and inspection is self explanatory. Cleaning 
involves the final cleaning of the finished assembly and any cleanup required during 
the assembly process, such as deburring of the holes after drilling. In some cases 
after drilling, it is necessary to disassemble the assembled parts for inspection and 
cleanup purposes. The disassembly process is essentially the reverse of the setup 
and clamping processes, and includes reassembly before the final securing is done. 

The program calculates a value in hours for each required subassembly and for the 
final assembly of the finished part. Currently, these calculations are based on 
constants stored with the BLOCK DATA. The actual equations take the following form. 

Setup Time: HOLES* TIME 6 Securing Time: HOLES * VALUE 

Clamping Time: CLAMP* TIME 1 Inspection Time: HOLES * TIME 2 

Drilling Time: HOLES* VALUE Disassembly Time: CLAMP* TIME 5 

Cleaning Time: HOLES * TIME 3 

where: HOLES = number of fasteners (fastener holes) 

CLAMP = number of clamps 

VALUE = volume of material removed from all the holes xa material 
complexity factor based on the material type xTIME 7 

The material complexity factor is assumed to be 3.0 for titanium. The constants 
TIMEX where X = 1 through 7 are stored in the BLOCK DATA as 

TIME1 = 0. 035 hour/clamp 
TIME2 = 0. 0008 hour /hole 
TIME3 = 0.01 hour/hole 
TIME4 = 0. 0015 hour /hole 
TIME5 = 0.012 hour/clamp 
TIME6 = 0. 0012 hour Aole 
TIME 7 = 0. 52 hour /cubic inch 

Provision has been built into the program so that the constants TIMEX can be replaced 
with equations. Values replacing the constants will be computed using these equations, 
which will be based on data collected and curve fit during the course of further study. 

An example of the type of data available is shown in Figure 2-16 where the hours 
required for drilling in titanium are shown as a function of the volume of material 
removed. 


2. 1.4. 4 Standard Hours . For each detail called out in the parts listing process, there 
is a corresponding list of required manufacturing operations in the cost subroutine. 

For each manufacturing operation a calculation of required standard hours is made. 
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Standard hours are defined as a standard 
time, measured in hours, representing 
an optimum required to perform a task. 
Standard hours for each production 
process are established by the industrial 
engineering department by analysis of 
time and motion studies of typical tasks. 
They are used as a means of measuring 
performance by determination of 
realization factors (or efficiencies) 
when compared with actual labor hours 
acquired through labor accounting 
processes. 


(cubic inches) 


Within the program, standard hours are 
calculated for two separate cost areas: 
factory production and final assembly. 
Factory production includes the manu- 
facturing and assembly of detailed parts 

into the basic subcomponents such as posts and beams. Final assembly involves the 
final assembly of the subcomponents into a finished TPS system. A listing of standard 
hours is output as part of the cost data under the heading STD HOURS (see Figure 2-13). 


Figure 2-16. Drilling Time in Titanium 


Standard-hour data have been collected for each of the seven manufacturing processes 
and assembly operations discussed; they have been adapted for use with this program. 
Standard-hour values are used to estimate the actual labor hours required for each of 
the production processes, and then to estimate the actual labor costs. Figure 2-17 
presents an example of a (Convair Aerospace San Diego operation) shop planning order. 
Listed are the various manufacturing processes required to produce a brace and the 
corresponding number of standard hours required for setup and running each shop 
process. The object of the standard-hour calculation technique in the program is, in 
effect, to be able to predict the planning order. 


Calculation of the standard hours for each of the manufacturing processes is performed 
in two parts: setup time and actual run time. Setup time is derived based on the size 
of the part, the complexity of the required setup, and the type of machine to be used. 
Run time is dependent on the machine operation rate, amount of material to be removed, 
depth of cut, and the surface area to be covered. The manufacturing standard hour 
equation is 


STDHR = SETUP + RUNTM * KT * SHIPSET 


where 

SETUP is the setup time in standard hours 
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GENERAL DYNAMICS CONVAIR ' " " ‘ _H^DAY~ AREA\_PA.RT NUMBER LT-N 

SHOP ORDER PLANNING 471 2 4T1 1015-133A A 
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06 
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001 

12 

050 
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001 

12 

060 

IDENTIFY TAG 
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065 

RW 













075 

INSP 
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05 
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ALK. ETCH STP57-301 
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ii 

090 

SULPH. ANOD STP5S-2C3 

5807 


0.0073 
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095 

INSPECT 




... - 







002 

03 

100 

C 11 ZC PRIMER STP59-201 

4815 
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0.0063 
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03 






110 

IDENTIFY R/S STP63-00]/ ' 
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0.07 

0.0032 
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TOTAL 


0.7 4 
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Figure 2-17. Example of a Shop Planning Order for a Brace 
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RUNTM is the run time in standard hours 

KT is the number of parts required of a given detail to produce a complete 
component 


SHIPSET is the number of shipsets of that component produced 


The standard data are usually derived for aluminum only. To apply the data to 
additional materials, material complexity factors are utilized. The material complexity 
factors account for the difference in manufacturing time requirements for performing 
identical tasks or operations on different materials. These factors are typically 
required only for those manufacturing operations associated with material removal 
such as drilling, milling, routing, burring, cutting, etc. Operations which usually 
do not require complexity factors are cleaning, layout, identification, painting, etc. 


Factors of 3.8 and 4.2 are commonly used as complexity factors for alloy steel and 
titanium, respectively. Methods are available for deviring complexity factors for 
other materials, but at present values based on the factors used in the COSTOT sub- 
routine are used. These factors are described under PROGRAM COST SUMMARY of 
this report. Table 2-3 summarizes the values for material complexity factors 
currently in use. 


Table 2-3. Summary of Material 
Complexity Factors Currently 
Being Used in the ASTDHR 
Subroutine 



Material 

Complexity 

Material 

Factor, F 

Aluminum 

1.0 

Titanium 

4.2 

Inconel 718 

4.2 

Hastelloy X 

4.0 

Rene 

4.2 

TD Nicrome 

5.0 

Columbium C-6752 

6.0 

Columbium C-129Y 

6.0 

Haynes 188 

4.2 

Tantalum 

6.0 

Alloy Steel 

3.8 

Stainless Steel 

4.2 

Magnesium 

4.0 

Beryllium 

5.0 

L-605 

4.2 


Calculation of the standard hours for each sub- 
assembly operation and for the final assembly is 
performed using the equations discussed in the 
previous section. The value of standard hours 
determined for subassembly processes are 
summed and added to the manufacturing process 
standard hours and listed for each detail part. 

The standard hours determined for the final 
assembly of the completed component are listed 
separately after each assembly process. 

2.1.4. 5 Realization . Realization factors are the 
ratio of standard hours assigned for a given process 
to the actual hours required. Realization for a 
process is used as a means of measuring actual 
performance against a standard time that is typically 
an ideal number of hours. It is a measure of shop 
efficiency and, as such, varies from department to 
department and from day to day within a department. 
Realization data for the various departments involved 
in production tasks at the Convair Aerospace San 
Diego operation have been collected, studied, and 
adapted for use with the program. Since these data 
take into account the effects of the learning curve 
for a specific task, realizations can be specified 
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either as an average value or as a time -dependent variable. Some of the factors 
affecting realization are 

a. Inaccurate planning of the required work, setup times, or run times. 

b. Machine breakdown. 

c. Change in machine, tools, or procedure which are not reflected by corresponding 
change in standard hour estimates. 

d. Tool breakage and part spoilage. 

e. Availability of previous setups. 

f. Use of special supervision. 

g. Ability and level of effort of individuals assigned to the task. 

Some of these factors are subject to control by managers and foremen, but they can also 
vary based on the current shop work load. Realization factors are useful in determining 
the overall effect of deviations from standard hour estimates. Realization for a 
particular task is defined as standard hours divided by actual labor hours. 

In the program, values for realization are stored for each manufacturing and assembly 
process, and since standard hours can be estimated, the actual labor hours can be 
calculated: 


Actual Labor Hours = Standard Hours/Realization 

Realization factors, along with labor rates and overhead ratios, are stored within the 
program. Provision has been made in the program for a different realization factor 
for each of the manufacturing processes, for subassembly, and for final assembly 
operations. However, at present, a constant realization factor of 0.40 is used for all 
operations, but data is readily available to establish an individual realization for each 
operation (although these factors would be strictly valid only within the Convair Aerospace 
San Diego operation shops). 

Average realization factors, along with the computed labor hours, are listed for each 
part in the cost output data under the headings REAL FACT, LABOR HOURS, and 
TOTAL HOURS. Figure 2-18 illustrates some typical realization factors and manufact- 
uring standard hours plotted as a function of the number of ship sets. 

2. 1.4. 6 Labor and Overhead Rates. Labor and overhead rates are used within the 
program to calculate appropriate costs, based on the number of actual labor hours 
required for each manufacturing and assembly process. 
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Figure 2-18. Typical Realization Factors and Standard Hours 

Labor rates reflect the wages paid directly to the individual employees for each hour 
of clock time. The rates do not include fringe benefits or company contributions to 
retirement, Social Security, and state unemployment; these are considered part of 
the overhead cost. Also included as part of overhead are indirect labor costs, 
maintenance, supplies, taxes, insurance, and depreciation. 

Labor rates are largely uncontrollable by management, being a function instead of 
union/management agreements and reflecting current labor supply and demand, 
general economic conditions, and inflation. Labor rates are a function of time and 
are readily predictable over the short term, although the incorporation of time- 
dependent equations has been left for future work. At present, the program can store 
a labor rate corresponding to each manufacturing process and the assembly operations. 
However, a constant manufacturing labor rate of $4.75 per hour is used, and for the 
assembly operations a constant rate of $5.20 per hour is used. 

The overhead ratio is the ratio of overhead costs to labor costs and can, therefore, 
be used to determine an effective overhead rate 

Overhead Rate = Overhead Ratio * Labor Rate 

The overhead ratio is a useful tool for estimating purposes, and is readily available 
based on past labor cost and overhead cost data. The program can store an overhead 
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ratio corresponding to each manufacturing process and the assembly operations. A 
constant ratio of 1.75 is used at present. Representative values of labor rates and 
overhead ratios as a function of time are presented in Figure 2-19. Values for the 
realization factor, labor rate, and overhead ratio are stored in the program. 

The average labor rate for each part is 
output with the cost data, along with the 
corresponding average overhead rate 
(overhead ratio for final assembly costs). 
Labor cost and overhead costs are calcu- 
lated for each part as the product of the 
labor hours and labor rate, respectively. 
Factory costs and assembly costs, then, 
are the sum of the corresponding labor 
and overhead coBts. 

2. 1.4. 7 Material Costs . Material costs 
are computed based on the material type 
and alloy (titanium, columbium, beryllium, 
etal), material form (sheet, plate, bar, 
et al) , and the raw material purchase 
weight. The actual calculation of material cost takes the form 

MATCOS = AMUV * KCOSWT * MAWT 

where MATCOS is the material cost in dollars for a given detail part, AMUV is the 
manufacturing usage variance factor explained below, KCOSWT is the material unit 
cost in dollars per pound, and MAWT is the raw material purchase weight described 
under WEIGHTS of this report. 

The calculation of material costs requires the derivation of a material unit cost KCOSWT 
and the definition of a material manufacturing usage variance factor AMUV (the calculation 
of MAWT is done during the weight analysis portion of the program). The parameters 
KCOSWT and AMUV originate in the AMATL subroutine, which is called from the ASTDHR 
subroutine once for each detail part, immediately prior to return to the COST subroutine. 
The actual material cost calculation takes place in the COST subroutine. 

The material unit cost KCOSWT is in general a function of the material type, alloy, and 
temper, material form, quantity of material bought, and special feature requirements 
such as special close tolerances, special lengths, widths, or thicknesses, and special 
marking. The determination of a value for material unit cost takes place in the AMATL 
subroutine where material price data has been curve fit and organized into equation form. 
Note that the integration of this subroutine into the program does away with the need for 
KCOSWT and KMUV K-TABLES previously provided in the BLOCK DATA subroutine. 



Figure 2-19. Typical Factory Direct Labor 
Rates and Overhead Ratios 
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Material form is specified by defining a value for KEY in the ASTDHR subroutine. 
Each detail part is assumed manufactured from the following forms: 


KEY - 0 

Insulation 

KEY = 4 

Wire, Rod, Car 

KEY = 8 

Casting 

1 

Fastener 

5 

Extrustion 

11 

Surface 

2 

Honeycomb 

6 

Tubing 


Insulation 

3 

Foil, Sheet, Plate 

7 

Forging 

12 

Ablative 


Material type is specified by input of a value for MAT LID , which represents the system 
primary structural material. The materials currently available in the program are 
summarized in Table 2-4. The defining of values for KEY and MATLID directs the 
program to a block of cost data within the AMATL subroutine which is associated with 
the given ma terial type and form. The material base price is then calculated utilizing 
the nominal material size requirements of length, width, and height (thickness) which 
are also carried over from ASTDHR with KEY and MATLID. 

Table 2-5 shows a portion of a typical price schedule for alloy steel plate between 0.25 
and 6. 00 inches thick. The equation resulting from a curve fit of the data takes the 
form 


Table 2-4. Summary of Primary 
Structural Material Types 
Currently Available in the Program 


Aluminum 

MATLID * 1 

Titanium 

2 

Inconel 718 

3 

Hastelloy X 

4 

Rene 41 

5 

TD Nicrome 

6 

Columbium Cb 752 

7 

Columbium C-1294 

8 

Haynes 188 

9 

Tantalum 

10 

Alloy Steel 

11 

Stainless Steel 

12 

Magnesium 

13 

Beryllium 

14 

L-605 

15 


Table 2-5. Example of a Typical 
Portion of a Material Price 
Schedule for Alloy Steel Plate 


Thickness 

E4340 AMS -6359 

(inches) 

Hot Rolled Annealed 

0.250 

45.00 

0.375 

44.40 

0.500 

44. 15 

0.625 

44.25 

0.750 

44.10 

1.000 

43.95 

1.250 

44.00 

1.500 

44.00 

1.750 

44.55 

2.000 

44.55 

2.250 

47.40 

2.500 

47.40 

2.750 

47.40 

3.000 

47.40 

3.500 

47.40 

4.000 

47.40 

4.500 

47.40 

6.000 

47.40 
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PBASE = . 006 * THK + .439 . 

where PBASE is the base price of the material in dollars per pound and THK is the 
thickness in inches. Thus, by specifying MATLID = 11, KEY = 3, and THK equal to 
the characteristic height, the program is directed to calculate the base price for a 
required size of alloy steel plate. In a similar manner the base price data for other 
material types and forms were curve fit and the resulting equations are available 
within the AMATL subroutine. 


The total material unit cost KCOSWT is equal to the material base price PBASE plus 
the sum of any extra cost item requirements PEXTRA. Table 2-6 presents a list of 
extra cost items available for aluminum plate, and Table 2-7 presents an example of 
one of the items, the quantity buy price differential. Typical extra cost items such 
as the quantity buy price differential, identification marking, packing for shipment, 
heat treatment, and ultrasonic tests, were curve fit in the same manner as the 
material base price data. They were organized in equation form for each material 
type and form to make up PEXTRA. In this way the total material unit cost + PEXTRA 
in dollars per pound is calculated directly within the AMATL subroutine. 


Table 2-6. Summary of Extra Cost 
Items Available for Aluminum 
Plate 


Actual Piece Count 
Alloys and Special Extras 
Circles 

Conversion Coatings 
Exact Quantity 

Identification Marking - Standard 
Identification Marking - Special 
Interleaving and Oiling 
Lengths, Long 
Lengths, Short 
Machined Surface (Two Sides) 
Mechanical Testing 
Packing 

Packing - per Mill-Std 649 
Sheet and Plate 
Protective Tape 
Quantity 

Test Material Samples 
Tolerances 

Diameter, Flatness, Length, 
Thickness and Width 
Ultrasonic Inspection 


Table 2-7. Example of the Quantity 
Buy Price Differential for 
Aluminum Plate 


Quantity per Item 
Pounds 

Extra 

Dollar /Pound 

300,000 and over 

Base 

29,999 - 

20,000 

0.010 

19,999 - 

10,000 

0.020 

9, 999 - 

8,000 

0.050 

7, 999 - 

4, 000 

0.070 

3, 999 - 

3,000 

0.125 

2, 999 - 

2,000 

0.285 

1, 999 - 

1,500 

0.450 

1,499 - 

1,000 

0.775 

999 - 

500 

1.775 

499 - 

250 

3.715 

249 - 

1000 

9.600 

99 - 

50 

19.800 
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For some combinations of material type and form, (e.g. , coated tantalum sheet) specific 
price data were not available. For these cases a characteristic material base price 
was established, as MBASE = 90.00 for coated tantalum. The specified material was 
then analyzed in terms of the equivalent aluminum material form (aluminum sheet), 
and the resulting value of PBASE derived for the equivalent aluminum form was ratioed 
using an aluminum base price (ALBASE = . 80) and the specified material base price 
(MBASE). Table 2-8 summarizes the values of MBASE in dollars per pound currently 
being used in the program. 

The manufacturing usage variance factor AMUV 
Table 2-8. Summary of Values for the is the ratio of the actual amount of material 
Characteristic Material Base Price purchased to the original estimated amount of 
Currently in Use in the Program material required for manufacturing. The 

factor is, in general, a function of the material 
form and past manufacturing experience. The 
factor results from material and part overbuy- 
ing to account for normal material losses 
during the manufacturing phase of production; 
it is established by accounting practice, as 
exemplified by Table 2-9. Material losses 
during manufacturing arise from part spoilage, 
loss, duplication, substitution, changes, and 
waste . 

A value for the AMUV is defined for each 
material form (for each value of KEY) in the 
AMATL subroutine. The defined value of 
AMUV and the calculated value of KCOSWT 
are returned to the COST subroutine for the 
actual calculation of material cost. At present 
a nominal or average value of 1.10, represent- 
ing a 10% overbuy, is used for all material 
forms. The final calculated value for 
KCOSWT is listed in the output for each detail part under MATL PER LB, and the 
resulting material cost under MATERIAL COST. The sum of FACTORY COST AND 
MATERIAL COST are listed as FABRICAT COST, and the total material costs for the 
complete system are summed under TOTAL MATERIAL COST. 

2. 1. 4. 8 Program Cost Summary. An example output of a program cost summary is 
presented in Figure 2-20. The total program costs are derived on the basis of the 
TPS gross weight, usually in terms of a Theoretical First Unit Cost or TFU. This is 
the predicted production cost of the prototype article including manufacturing plan n ing, 
fabrication, subassembly, sustaining engineering, sustaining tooling, quality control, 
materials, and subcontracted parts. 
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Table 2-9. Typical Manufacturing Usage Variance Factors for a 
Past Commercial Transport Aircraft Program 


Contract 
Lot No. 

Material 
Costs, A 
(millions of $) 

Original Estimated 
Material Costs, E 
(millions of $) 

Percent Variance 
(A - E/A) * 100 
(percent) 

Mfg. Usage 
Variance 
Factor 

1 

40.65 

34. 74 

17.0 

1.170 

2 

4.61 

4.25 

8.5 

1.085 

3 

16.67 

14.39 

13.8 

1.138 

4 

22.69 

21.40 

6.1 

1.061 

5 

16.28 

15.84 

2.8 

1.028 

6 

66.50 

62.15 

7.0 

1.070 

7 

10.22 

9.84 

3.9 

1.039 

8 

68.71 

61.94 

10.9 

1.109 


THERMAL PROTECTION SYSTEM COST SUMMARY 


COST $/3Q FT 

THEORETICAL FIRST UNIT COST (TFU) 461.49 


NON RECURRING COST 


ED AND D 

22248.04 

TOOLING 

1172.26 

GROUND TEST HARDWARE 

1730.58 

FLIGHT TEST ARTICLES 

922.98 

FLIGHT TEST S AND RP 

309.20 

TOTAL NONRECURRING TPS COST 

26383.05 

RECURRING PRODUCTION COST 


SUSTAINING ENGINEERING - INCLUDED IN TFU 
SUSTAINING TOOLING - INCLUDED IN TFU 

PRODUCTION ARTICLES 

461.49 

TEST ARTICLE CONVERSION 

138.45 

TOTAL RECURRING PRODUCTION COST 

599.93 

RECURRING OPERATIONS COST 


REPLENISHMENT S AND RP 

614.70 

TOTAL RECURRING OPERATIONS COST 

614.70 

TOTAL TPS PROGRAM COSTS 

27597.68 

NUMBER OF PRODUCTION UNITS 

1.0 


Figure 2-20. Program Cost Summary 
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The following expressions develop thermal protection system total cost. This is a 
weight-driven procedure which employs a technique independent of those previously 
discussed. Values of the embedded constants are given in Table 2-10 and complexity 
factors are given in Table 2-11; however, specification of the program cost option 
ICOS as unity allows the user to input his values for these parameters. 


Theoretical first unit cost (TFU) 

TFU = CF X K x (W) Ai 

Engineering Design and Developments Costs (EDD) 
EDD = K 2 (W) bt 
T ooling Costs (TOOL) 

TOOL = K 3 (W) CT 
Hardware requirements 

Ground test hardware (GTH) 

GTH = EGTH TFU 
Flight test articles (FTA) 

FTA = FTH TFU 


| These parameters 
r determine total 
1 non-recurring costs. 


Flight test spares and replenishment parts (FTSRP) I 
FTSRP = FTS TFU J 


Production articles(PA) 

PA = NPA TFU 
Test article conversion (TAC) 

TAC = ETA TFU 

Replenishment spares and replacement parts 
RSRP = OS * TFU 


■\ These parameters, along 
with sustaining engineering 
„ and tooling (both included in 
the TFU) determine total 
recurring production cost. 

* 

(RSRP) 1 This parameter determines 
V total recurring operation 
I cost. 


The total program cost (TPC) is obtained by a summation of non-recurring DDT&E(NR), 
recurring production (RP), and recurring operations (RO) costs. 

THEORETICAL FIRST UNIT COST. The parameter which is fundamental to computing 
all other terms of the total program cost is the theoretical first unit cost (TFU). This 
is the predicted production cost of the prototype article including manufacturing, 
planning, fabrication, subassembly, sustaining engineering, sustaining tooling, quality 
control, materials, and subcontracted parts. 
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Table 2-10. Values of Table 2-11. Summary of the Complexity Factors Currently 
Total Cost Constants Being Used in the COSTOT Subroutine 


AT 

0.667 


Material 

Structural 








Type 

Configuration 

Complexity Factor 


BT 

0. 187 


MATLID 

KINDP 


CF 



CT 

0.86 


1 

1 2 3 4 5 1.3 

2.1 

1.0 

2.1 

2.1 

K1 

0.00171 


2 

3.5 

6.0 

3.0 

6.0 

6.0 




3 

3.2 

4.3 

2.5 

4.3 

4.3 

K2 

6. 58 


4 

3.1 

3.1 

1.9 

3.1 

3.1 

K3 

0.004324 


5 

4.0 

4.7 

3.0 

4.7 

4.7 




6 

5.2 

6.9 

4.2 

6.9 

6.9 

EGTH 

3.75 


7 

17.0 

17.0 

12.0 

17.0 

17.0 

FTH 

2.00 


8 

17.0 

17.0 

12.0 

17.0 

17.0 




9 

3.2 

4.3 

2.5 

4.3 

4.3 

FTS 

0.67 


10 

19.0 

19.0 

14.0 

19.0 

19.0 

NPA 

1.00 


11 

1.8 

2.6 

1.4 

2.6 

2.6 




12 

2.2 

3. 0 

1.7 

3.0 

3.0 

ETA 

0.30 


13 

2.2 

3.0 

1.7 

3.0 

3.0 

OS 

1.332 


14 

6.0 

10.0 

5.0 

10.0 

10. 0 ; 




15 

3.8 

4. 6 

2.8 

4.6 

4.6 i 

STPS 

22,000 









NONRECURRING DDT&E COST. This program cost consists of the engineering design 
and development (ED&D) cost and the production cost for ground test hardware (GTH), 
flight test hardware (FTH), and flight test program spares (FTS). The ED&D cost 
includes the engineering, design, development, laboratory test, support activities for 
subsystem development, and the cost of hardware required to support laboratory 
development testing and component qualification testing. Initial tooling for the fabrica- 
tion of the subsystem is included in ED&D. Ground test hardware costs reflect all 
hardware used in both subsystem and combined subsystem development testing. Hard- 
ware costs are generated by multiplying the equivalent number of units required by the 
theoretical first unit (TFU) cost of the subsystem. No cost improvement (learning) is 
taken during the development program. The TFU cost of the subsystem includes 
manufacturing planning, fabrication, subassembly, sustaining engineering, sustaining 
tooling, quality assurance and control, and materials and subcontracted parts. 

RECURRING PRODUCTION COST. This program cost covers the cost of manufacturing 
the production articles (PA) and converting the test articles (TAC) to the operational 
configuration. These costs are generated by multiplying the equivalent numbers of units 
by the TFU. No cost improvement (learning) is taken. 

RECURRING OPERATIONS COST. This program cost covers the cost of spare hard- 
ware required during ten years of operation. The cost is generated by multiplying the 
equivalent subsystem replacement fraction for 444 flights (OS) by the TFU. Again, no 
cost improvement (learning) is taken. 


2-36 





The derivation and supporting rationale of the costs for basic subdivisions of work, i.e., 
engineering (ED&D), manufacturing (TFU and hardware requirements), and tooling, 
are generally based on one or all of the following: (1) historic cost data actuals, (2) 
vendor cost data, with or without adjustments, received in support of this or other 
studies, and (3) other cost models. 

Since these nhodels are all based upon a total TPS weight, all computations are made 
for a total weight based upon a total area of the order of 10,000 ft^. Hence, total 
weight is simply the product of panel unit weight and system area, and panel unit cost 
is derived from the quotient of system total cost and system area, the dividend driven 
by system total weight. The number of hardware units, flight test spares, and the 
like are then really portions of the total vehicle TPS and not simply an additional 
number of panels. 

2. 1.4. 9 Total Program Cost Analysis for RSI, Ablators, and Carbon- Carbon . An 
additional task was to establish a method for assessing the relative cost differences 
between alternate, orbiter-vehicle thermal protection systems (TPS). A parametric 
cost model was developed which calculates TPS total program costs in three 
categories: nonrecurring (DDT &E), recurring production, and recurring operations 
(maintenance and refurbishment). This concentrated on three specific TPS concepts, 
which were reusable carbon- carbon composites (C/C), reusable surface insulation 
(RSI), and ablatives. 

Approach 

Since the primary application of costs in this program is for design selection and 
optimization, emphasis is placed on the determination of relative cost differences 
rather than absolute costs such as those normally used for contract bidding purposes. 
The approach used in developing the parametric cost model started with the collecting 
of available cost data on thermal protection systems. A work breakdown structure 
(WBS) of major elements of cost was then established for typical thermal protection 
systems. Next, a set of parametric cost estimating relationships (CERs) were 
developed which correlated costs with weight for the various WBS elements. A sample 
calculation demonstrating the use of the model was made for a hypothetical orbiter TPS. 

The cost data that were collected for the analysis could be grouped into three general 
categories, (1) actual costs from historical programs, (2) estimated costs from 
proposals and study results, and (3) existing cost models. 

The actual costs from historical programs were of course the most desirable data 
because they are most representative of true program costs. However, these data 
were mostly metallic concepts and consequently only of limited use as a source of 
analogs for the C/C, RSI and ablative concepts being analyzed in this study. The 
costs appearing in recent space shuttle proposals and study reports were more similar 
to the TPS designs of interest in this study, but were estimates and not actual costs. 
Data from these sources typically suffer from a lack of design definition (e.g. .uncertain 
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weights, no m aterial distributions, no areas, locations, etc), inadequate cost visibility 
(e.g. , no unit manufacturing costs shown, learning assumptions unknown, TPS costs 
buried in structure, etc), or general incompatibility between the design data shown and 
the cost data. Existing cost models used by Convair and other contractors were 
analyzed and utilized in the study. The actual data points which these models use are 
usually not shown due to proprietary reasons, but general trends such as the slope of 
various CERs and knowledge of a contractor's product line provided some data for our 
analysis. Some of the data sources utilized in this analysis were 

General Dynamics Historical Data 
General Dynamics Space Shuttle Cost Model 
North American Rockwell Space Shuttle Cost Model 
North American Rockwell Historical Data 
Aerospace Corporation Earth Orbit Shuttle Cost Model 
Space Shuttle Phase B and Phase B" 

McDonnell Douglas Corporation Cost Model 

The work breakdown structure established for the cost model is shown below. 

TPS Total Program Cost xxx 

Nonrecurring xx 


Development 

XX 


Test Hardware 

XX 


Tooling 

XX 


Recurring Production 


XX 

Production and Test Vehicle Mods 

XX 


Recurring Operations 


XX 

Operational Spare and Repair Parts 

XX 



The elements of cost shown were based on an analysis of the space shuttle final reports 
and are representative of the total cost for a typical orbiter vehicle thermal protection 
system. 

The development of the parametric cost estimating relationships (CERs) was the major 
task in building the cost model. Theoretical first unit (TFU) production cost CERs were 
generated for RSI, C/C, and ablative type TPS concepts. CERs for typical attachment 
structure and for insulation material were also developed because these two components 
were generally required in most TPS designs. A similar set of CERs for engineering 
design and development was also developed. A CER for initial tooling costs was 
developed as a function of vehicle total TPS weight and was independent of the type of 
TPS being evaluated. The refurbishment hardware costs for recurring operations 
were developed in the form of a range of cost-per-flight factors for each of the three 
TPS concepts. The range in these factors for any given type of TPS reflected the 
specific application (severity of environment) or location on the orbiter vehicle. 
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In the development of CERs, statistical regression techniques would normally be applied 
to a scatter plot of historical data points; however, in this study the lack of sufficient 
historical data prohibited the use of this method. A modification of this technique was 
used in which the slope of the CER equation was determined from other supporting data 
or existing CERs; then a line with this slope was passed through the best available 
analogs for the component in question. In other cases, CERs were adopted from 
previously developed cost models if they appeared to be in good agreement with the 
data and with similar CERs from other sources. 

The sample calculation utilized a combination of fully reusable TPS types. Program 
characteristics (numbers of inventory vehicles, number of flights, test hardware 
conversion assumptions, etc) used in the sample calculation were based on the Space 
Shuttle program Phase B study period. 

The followi ng ground rules and assumptions were utilized in the analysis. 

1. The development cost includes the engineering, design, development, laboratory 
test, support activities, and the cost of hardware required to support laboratory 
development testing and component qualification testing. 

2. The development hardware costs represent the production cost for ground test hard 
ware, flight test hardware, and flight test program spares. 

3. The tooling cost includes all activities for engineering and building of initial tooling 
required for TPS fabrication. 

4. The recurring production cost includes the cost of manufacturing the production 
articles, converting test articles to operational configuration, and sustaining 
tooling activities. 

5. The recurring operations cost covers the cost of spare and repair hardware 
required for the flight program. 

6. An operational program of 445 flights was assumed. 

7. The following hardware requirements expressed in terms of equivalent TFU costs, 
were assumed: 

Development test hardware consumed = 1. 7 
Production articles and flight test conversion = 3.6 

8. All costs are in terms of 1970 dollars and exclude fee. 

9. No costs are included for vehicle-level activities such as flight test operations, 
ground support equipment procurement and maintenance, etc. 
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Results 


The cost model results are summarized in equation form in Table 2-12. In the following 
paragraphs the development of each CER is discussed. The reader should remember 
that the parametric costing approach required in a study of this type is by nature depend- 
ent on historical cost experience. The quality and quantity of historical data on thermal 
protection systems analogous to the Space Shuttle orbiter vehicle are very limited. 
Additionally, the relationship between TPS cost and TPS performance is not well under- 
stood at this time; consequently, in many cases a single previous data point or trends 
observed from historical airframe data were utilized in generating CERs. 

TFU CERs 

The theoretical first unit (TFU) manufacturing cost CERs (Figure 2-21) were based on 
an analysis of the North American Phase B study results, and existing cost models 
from Convair Aerospace, Aerospace Corporation, and McDonnell Douglas. A cost/ 
weight scaling exponent of 0. 6 was selected from the Aerospace model. This value is 
in good agreement with historical total airframe data that show cost/weight scaling 
exponents ranging from 0. 5 to 0.8. The McDonnell Douglas model CERs are a function 
of area and individual panel size; therefore, no direct comparison of selected scaling 
exponent could be made. The RSI and carbon-carbon concepts are shown as having 
the same range of costs depending on the specific application. The upper edge of the 
band shown in the CER would correspond to a more complex application such as the 
wing leading edge and nose cap areas. The lower edge represents a general application 
over more gently contoured and lower temperature areas on the remaining body areas 
and aerosurfaces . 

No good evidence was found for distinguishing between the two concepts, and since 
both are still in the development stages, any claims of "lower cost" must be considered 
suspect. Reported costs from individual proponents are often based on optimistic or 
exclusive applications that may not be realized in initial operational configurations. It 
seems likely that a significant portion of each type will be used on the first orbiter 
vehicles with carbon-carbon being used in the high-temperature highly contoured areas, 
RSI in the low-temperature gently contoured areas, and both competing for the 
moderate-temperature, moderately contoured areas. 

The ablative CER is shown at about one-half the complexity of the RSI and C/C. This 
line passes half way between two independent estimates of an ablative TPS design, and 
the ratio of complexities between the ablative and advanced (RSI & C/C) concepts is in 
rough agreement on a similar cost model . 

The attachment structure CER was based on the Aerospace cost model and assumed the 
use of titanium for this application. The North American orbiter TPS carrier panels 
estimate fell just above (more expensive) this parametric line as could be expected due 
to the use of a Ti sandwich construction. The insulation CER was taken directly from 
the Aerospace cost model. 
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Table 2-12. TPS Cost Model Summary 
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♦RSI and C/C TFU la represented by a band rather than a single line. 
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Initial Tooling CER 

The initial tooling cost CER was selected from an existing set of tooling CERs developed 
for the Space Shuttle booster vehicle (Figure 2-22). The particular line selected was 
one used for estimating aerosurfaces initial tooling requirements. The physical shapes 
and contours encountered in aerodynamic surfaces were felt to closely approximate the 
TPS shapes and contours. Values obtained from this CER were in good agreement with 
those from the North American final report, McDonnel Douglas thermal-structure 
tooling CER, and the Aerospace cost model. 

ED&D CERs 

The engineering design and development CERs (Figure 2-23) used the slope from the 
Aerospace model TPS. A line with this slope was then passed through the North 
American orbiter development estimate for an RSI type TPS with C/C nose and wing 
leading edges. No difference between RSI and C/C was indicated and this CER is 
properly entered at the total weight of the RSI and C/C. The attachment structure 
CER was established by passing another line of the same slope through the estimated 
development cost for the titanium carrier panels on the NAR orbiter. This line 
coincided with the Aerospace cost model's baseline (complexity factor of 1. 0) thermal 
protection panel CER. The ablative TPS development CER was then estimated to have 
a relative complexity of about 1.2 based on the Aerospace model's allowance for 
materials and configuration complexities. This resulted in a line that was about 50% 
of the RSI and C/C line as was the case for TFU. The CER for TPS insulation was 
taken directly from the Aerospace model even though this was about 25% higher than 
the NAR orbiter insulation data point. 

Refurbishment CERs 


The two contributors to TPS refurbishment costs are hardware consumption, in the form 
of spare and repair parts, and labor. For this analysis the differences in labor require- 
ments between the various TPS concepts were assumed to be negligible and were excluded. 
This assumption is reasonable considering the objective of obtaining relative cost 
differences rather than absolute costs for a specific concept. Historically, the labor 
portion of operations cost is small compared to the materials used. One reason for 
this is that labor is not so "dedicated” as hardware. A mechanic can work on landing 
gears, wings, etc when he is not doing TPS work (i.e. , more efficient utilization of a 
dollar spent on labor versus hardware). Also, this labor portion is unlikely to be 
greatly different between TPS types because a large part of the labor task is associated 
with inspection to determine which areas are qualified for an additional flight, and this 
must be done regardless of the design concept. Removal and replacement activities 
will be relatively insensitive to TPS type. Here the more important discriminator is 
likely to be the panel size with smaller panels requiring more labor due to the greater 
number of dissimilar interfaces which produce problems with joint mismatch, seals 
and plugs, etc. The scope of this study did not permit going into these details for each 
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Figure 2-22. Initial Tooling Costs 
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type of TPS. The hardware costs are simply a function of the quantities consumed and 
their respective unit costs. This consumption was expressed in terms of equivalent 
theoretical first units per flight, which are multiplied by the TFU costs to get operational 
spares costs. As in the other areas of the cost analysis, the historical data base 
applicable to orbiter TPS consumption is minimal. Additionally, the replacement 
rate of any TPS concept will likely be affected by the degree to which reusability is 
stressed during development. We can also expect the replacement rate to change 
(either increase or decrease) through the life of the program. For example, removal 
frequencies of certain panels may be changed after the initial flights because the 
predicted thermal loads differ from actual experience. The specific application would 
also affect the consumption rate. For example, ablatives used on the nose and wing 
leading edge may require 100% replacement per flight, while on other areas of the wing 
and body they may withstand 10 or 20 flights. The range of spares consumption rates 
considered for the TPS concepts analyzed is shown in Figure 2-24. 

Sample Calculation 

The calculation of a typical TPS total program cost is shown in Figure 2-25. This 
TPS concept is similar to the design submitted for the North American Phase B 161-C 
orbiter and assumes the same hardware equivalents for the program. 

Uncertainties 


The results of the economic analysis must be considered uncertain because of a number 
of factors peculiar to this study. 

1. No schedule effects were considered; consequently, time phasing of costs is not included 
in the selection criteria. Benefits associated with lower peak-year funding requirements 
or smaller early-year development outlays are not defined by the cost model . In 

view of the current 
budget constraints in 
aerospace and related 
fields, this is a limita- 
tion of the study. 


2. The general absence of 
actual historical experi- 
ence in developing, 
manufacturing, and 
operating systems such 
as those investigated in 
this study makes the cost 
model results of this 
analysis highly uncertain. 
Adding to this uncertainty 
Figure 2-24. TPS Spares Consumption Ranges (mainly in the CERs 
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developed) is the fact that the designs being evaluated are still preliminary in 
nature with relatively little detail. One would expect a significant difference 
between the description of any TPS at this time and the actual operational system 
that results. The cost model, however, is intended to predict the cost of the 
operational system without the benefits of knowing what the final design will be. 

3. Existing cost models have not dealt with these types of systems except on a purely 
theoretical basis (i.e. , "paper" designs). Differences inherent in the existing cost 
models, such as different cost categories, different cost driving parameters, 

(e.g. .weight vs area), and different ground rules (no fee, no sustaining engineering, 
no program management, etc) prohibited making a comparison of results on a 
single TPS concept. Since this analysis relied heavily on current cost model 
experience, the validity of the results must be considered uncertain. 

Values of the total program cost parameters have been developed here (Table 2-13) to 
allow the program user to input applicable parameters. 

2. 1. 5 MATERIAL PROPERTY DATA . A limited effort was expended throughout the 
contract to gather and reference the material property data necessary to run the TPS 
sizing procedure. Much of the data gathered in the previous contract, NAS9-10956, is 
of course applicable to the metallic cover panel for which the majority of that effort 
was spent. The efforts of the present contract have been directed to four new systems - 
RSI, ablators, carbon-carbon, and phase change materials. The following paragraphs 
discuss references and what properties are available for each material. 

The properties of metallic s 
have been amplified by the 
publication of data by General 
Electric's publication, Space 
Shuttle Thermal Protection 
System Metallic Materials 
Property Data (Ref . 8). The 
metallics considered are 
outlined in Table 2-14 along 
with the properties included 
in the handbook. It should be 
noted that not all properties 
are necessarily available, 
but at least an attempt is 
made to provide them. 

Unfortunately, no single 
compendium exists for 
properties of the other three 
concepts: RSI, ablators, and 


Table 2-13. Total Cost Program 


Param- 

eter 

NS] 

ECT 

1-6 

7 

8 

9 


RSI 

Ablator 

C-C 

CF1 

Table 

1.0 

1.0 

1.0 

AT 

0.667 

0.6 

0.6 

0.6 

BT 

0.187 

0.6 

0.6 

0.6 

CT 

0.86 

0.86 

0.86 

0.86 

K1 

0.00171 

0.06528 

0.03295 

0.07181 

K2 

6.58 

0.575 

0.300 

0.575 

K3 

0.004324 

0.004974 

0.004974 

0.004974 

EGTH 

3.75 

1.7 

1.7 

1.7 

FTH 

2.0 




FTS 

0.67 




NPA 

1.0 




ETA 

0.3 




OS 

1.332 




STPS 

1000 

1000 

1000 

1000 
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Table 2-14. List of Metallic TPS Property Data Contained in Reference 8 


Materials 

Data 

Inconel 718 

Commercial Designation: 

Mechanical Properties: 

Ronp 41 


Elastic Modulus 

TD-Nickel Chromium 

Composition: Weight Percent 

Shear Modulus (or 

(2% Th02) 

Physical Properties: 

Poisson's Ratio) 

L-605 (Haynes Alloy 

Thermal Expansion 

0. 2% Yield Strength 

No. 25) 

Thermal Conductivity 

Ultimate Tensile Strength 

Haynes Alloy No. 188 

Specific Heat 

Tensile Ductility 

Cb-752 

Emittance 

Creep 

FS-85 

Density (Room Temp) 

Creep Rupture 

C-129Y 


Bend Ductility 

B-66 


Impact Properties 

T-lll 


High -Cycle Fatigue 

T-222 


Low-Cycle Fatigue 


carbon-carbon. The reusable surface insulation system has been studied in great 
detail by three companies: General Electric, Lockheed, and McDonnell Douglas. 

The final reports of these companies' efforts contain all the current physical and 
mechanical property values. Since the reports themselves are available in either 
hard cover or microfiche, only the references are given here (References 9, 10 and 11). 
Both test and design values are given for all RSI properties as well as data for the glue 
line, the substrate, and other related parts of the configuration. It was necessary to 
select a consistent set of RSI properties to be used in calculations during the course 
of the study. For this, the General Electric data of Reference 11 was selected as 
typical, and it appears in Table 2-15. This selection was based primarily on the ease 
of computer implementation considerations and does not imply a preference of 
recommendation for this material over those supplied by other contractors. 

Typical property values for ablative materials are usually available from the material 
manufacturers. Since no structural analysis is performed on the ablative system, only 
thermodynamic properties are of interest here. Table 2-16 gives the thermal, 
thermophysical, and mechanical /physical properties of two charring ablators 
manufactured by General Electric. These are typical, but if other property data is 
desired, the program user is urged to check with ablator manufacturers. 

Material property data for carbon-carbon composites are difficult to obtain basically 
because of limited industrial experience in dealing with the materials. The major 
efforts to date have been performed by the Vought Missiles and Space Company in 
their Phase II contract on the TPS development for the Space Shuttle wing. Other 
companies, including Convair Aerospace, have been developing carbon-carbon concepts 
under company funded research or contracts to other agencies such as NRL. To 
obtain data for the stringer reinforced leading edge configuration, Reference 12 is 
recommended. For the truss core configuration, the thermal and mechanical properties 
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Table 2-15. Typical Properties of Reusable Surface Insulation (RSI) 


TEMPERATURE 

THERMAL CONDUCTIVITY 
(Btu/ft sec R) x 105 

SPECIFIC HEAT 
(Btu/lbm R) 

R 

F 

P < 

10" 4 ATM 

P = 

lO- 2 ATM 

P = 

1 ATM 

REI MULLITE 


40 

0.187 

0.550 


0.172 



0.215 

0.610 

ar-aaw 

0.195 

700 

1 " 

0.248 

0.673 


0.214 

800 


0.283 

0.743 

1.20 

0.228 

900 

440 

0.33 

0.820 

1.30 

0.240 


540 

0.376 

0.910 

1.41 

0.250 


790 

0.532 

1.16 

1.72 

0.267 


1040 

0.756 

1.49 

2.1 

0.280 


1290 

1.07 

1.90 

2.57 

0.289 

1 

1540 

1.60 

2.45 

3.14 

0.295 

2250 

1790 

2.14 

3.13 

3.86 

0.300 

2500 


3.03 

4.00 

4.71 

0.300 

2750 


4.30 

5.15 

5.80 

0.300 

3000 

2540 

5.10 

6.60 

7.07 

0.300 

DENSITY (lbm/ft 3 ) 

15.0 



developed at Convair Aerospace under company funded research are presented in 
Tables 2-17, 2-18, and Figure 2-26. The program user is cautioned however, that 
material properties are extremely process and configuration dependent - to ensure 
their proper use, only property values ascertained by test for the particular configura- 
tion should be used. 

For the case of the heat sink TPS (the phase change material - PCM), used in 
conjunction with the RSI system, thermodynamic properties for a wide variety of 
PCMs have been gathered at Lockheed Huntsville (Reference 13). These fillers 
include both high and low temperature phase change materials as well as a number 
of solid-solid state change ones too. Additional information on heat sink systems is 
being developed by McDonnell Douglas under Contract NAS8-27708. A limited 
amount of data has also been gathered to describe the thermodynamic properties 
of several fiberglass and aluminum honeycomb configurations. These concepts are 
applicable to the honeycomb substrate which will be used for the RSI both with and 
without an impregnated PCM. The honeycomb properties are equivalent specific 
heat and thermal conductivity, the latter being evaluated at low temperature at which 
radiation effects across the substrate are negligible. These data will be applicable to 
the RSI TPS since the thermal constraint on the system is a bond line temperature on 
the order of 300F . The thermodynamic properties are presented in Figures 2-27 
through 2-32, and Tables 2-19 through 2-23. 
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Table 2-16. Properties of Low-Density Elastomeric Shield Materials (ESM 1004) 


PROPERTIES 


ESM 1004-X 


ESM 1004 AP 


THERMAL 

A. Specific Heat (Btu/lb-F) 

OF 
100F 
3 OOF 
500F 

B. 

Dynamic Enthalpy 


C. 

Thermal Conductivity x 10 5 

(Btu/ft-sec-F) 


1 atm, air 

OF 



100F 



300F 



Strong Endo therm ■ 
0.01 atm, air 
Slight Exotherm - 
1 atm, air 


Exotherm 

(Air, 1 atm, 725F) 



D. TGA - 5% Weight Loss *650F 800F 

50% Weight Loss 1150F 1120F 

Residual Weight 35% 43% 

*ESM 1004A 0 Post Cure = 500F, ESM 1004X Post Cure = 300F. Initial Weight Loss Due to Unreacted 
Residual Volatiles and Not Thermal Decomposition \ 


E. Weight Loss in Vacuum 

10' 5 MM HG - R.T. 5 Days 


THERMOPHYSICAL 

Virgin Density, p v lb/ft^ 

Char Density, P c , lb/ft® 

Surface Emissitity 

Pyrolysis Gas Specific Heat, CPg Btu/lb R 
Molecular Weight of Injected Species, Mq 
O rder of Reaction 
Pre-Exponential Factor, Z sec -1 
Activation Energy, E Btu/lb mole 
Heat of Decomposition Hqf 1135R 

Btu/lb Gas Generated 1460 

1710 

2610 


MECHANICAL/PHYSICAL 
A. Density (lb/ft 3 ) 



B. Tensile Properties 





Elastic 




Stress 

Elong. 

Modulus 

Stress 

Elong. 


(psi) 

(%> 

(psi) 

(psi) 

(%) 

-150F 

32 

40 

95 

150 

50 

-100F 

14 


— 

90 

46 

-50F 

9 

22 

60 

70 

34 

75F 

7 

22 

45 

60 

25 

300F 

6 

17 

45 

50(est) 

22 


Elastic 

Modulus 

(psi) 



C . Bond Shear Strength 
Shield/Bond System 


Shear Strength (psf) 
(RTV 560 to Ti) 



D, Bond Shear Strength - Bond Aluminum to 
Aluminum RTV 560 Adhesive 

75F 

400F 


Thermo Structural 
Compatibility 


Low Temperature 


Shear Strength (psi) 


600 (Average) 
250 (Average) 


No Failure Down to -300F 
as Substantiated by Test 


Shear Strength (psf) 
(RTV 560 to BE) 


(Cohesive Failure 
in ESM 1004AP) 


Shear Strength (psi) 


600 (Average) 

250 (Average) 


No Failure Down to -300F 
as Substantiated by Test 


































Table 2-17. Assumed Material Properties 


Carbon- Carbon (Pyrocarb 500) 

Density = 93. 7 lb x ft" 3 

Specific Heat =0.22 Btu x lb - * x (R) -1 

Thermal Conductivity* in Btu x ft"! x hr"! x (R)“! 


Temperature, R 

960 1460 1960 2460 

Conductivity 

.958 .992 1.03 1.06 


Dynaflex Insulation 

Density = 1 lb x ft -3 
Specific Heat in Btu x lb"! x 


Temperature, 

R 

500 

1000 

1500 

2000 

2500 

3000 

Specific Heat 

.177 

.220 

.242 

.257 

.270 

.277 

Thermal Conductivity in Btu x ft"! x hr"* x 

(R)" 1 



Temperature, 

R 

660 

1060 

1460 

1860 

2260 

2460 

Conductivity 

.0225 

.0417 

.0625 

.0834 

.121 

.141 


Titanium Structure 

Density = 281 lb x ft" 3 

Specific Heat = . 14 Btu x lb” 1 x (R)"! 

Thermal Conductivity = 5.0 Btu x ft"! x hr -1 x (R)" 1 
♦Taken same as for Pyrocarb 400 


2-52 




















2-53 

























lo o ld © m © 


• ••••• 

CQ (M rH iH © © 

x t _hh x T _xa x nxa 
AXIAIIOnaNOO 1V]/\EH3HX XN UHVddV 


2-54 


Core 








0.4 



h ‘qi/nia ‘xvaH Diaioaas 


2-55 








0.12 



h g -tVnxa ‘AxiAixonaNoo ivkh^hx 


2-56 


Figure 2-28. Thermal Conductivity of Fiberglass Honeycomb Sandwich Panels 
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Figure 2-30. Specific Heat of Honeycomb Sandwich Panels Having Alu minum Core 
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Figure 2-32. Thermal Conductivity of Aluminum Honeycomb Sandwich 


Table 2-19. Specific Heat of Fiberglass Honeycomb 
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2.1.6 AEROTHERMODYNAMIC ANALYSIS . The prediction of aerothermodynamic 
environment is based on the Convair Aerospace Computer Program P3020 (Ref. 14 
and 15) and the recent recommendations of the NASA thermal panel on heat transfer 
methodology (Ref. 16). Figure 2-33 outlines the computational procedure. Using the 
Patrick AFB atmosphere, the free stream conditions are determined for a given flight 
trajectory. The boundary layer edge conditions are then established by the prediction 
of the flow field about simple geometric shapes such as flat plates, wedges, cones, 
cylinders, or spheres. Cross-flow theory is then applied to correct the three- 
dimensional effects. 

Prediction of flow field is classified into two regimes: high and low local angle of 
attack. For low angle of attack, the shock waves are assumed attached to the body, 
and flow properties are computed from tangent wedge/cone techniques. At high angles 
of attack, the flow field cannot be predicted so conveniently as at low angles of attack 
because of shock detachment. Thus, approximate methods are used. 

Computation of local heating rates is classified into three regimes: laminar, turbulent 
and transitional boundary layer. For the laminar boundary, the Blasius solution iB 
evaluated at the Eckert reference enthalpy, and heat transfer is computed from Colburn's 
statement of Reynolds analogy. For turbulent flow, two techniques are presented; the 
first utilizes skin friction predicted by the Schultz -Grunow equation evaluated at Eckert's 
reference enthalpy, with Colburn's Reynolds analogy, and the second is the Spalding-Chi 
skin friction and von Karman's Reynolds analogy. Masek's criterion is used to predict 
the onset of transition. Transitional heating between the laminar and turbulent boundary 
layers is then calculated as a linear interpolation of turbulent and laminar heating values, 
the degree of turbulence depending on the "turbulent fraction" exhibited by the boundary 
layer with respect to values of Reynolds number for transition onset and end. The 
virtual origin for turbulent flow can be taken as either the leading edge or the beginning 
of transition. Swept cylinder theories of Fay-Riddell (laminar) and Beckwith-Gallagher 
(turbulent) are used to compute heating rates on leading edges and for bodies at high 
angle of attack. 

Flow Field Calculations 

The effective angle of attack, is defined as the scalar product of the direction cosines 
of the normal to the surface and of the velocity vector. Using Figure 2-34 these are 

x, y, z direction cosines of normal to surface 

DNX = sin 6, DNY = - cos 9 cos 0, DNZ = - cos 9 sin 0 

x, y, z direction cosines of velocity vector 
Vx = V cos a cos /3, V z = V sin a cos j9 
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Aerodynamic 

Heating Prediction Transition Criteria Flow Field Properties 



IC0NF = 1 Flat Plate With Cross Flow 

2 Cylinders 

3 Cone 

4 Wedge 

5 Sphere 

6 Laminar Cylinder 

7 Turbulent Cylinder 


Figure 2-33. Aerodynamic Heating Prediction 
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(2-23) 


Then, the scalar product is 


sin f = DNX cos a cos j8 + DNY sin jS + DNZ cos 0 sin a 


where 

a- = angle of attack 
jg = yaw angle 

Thus, body orientation can be described 
simply as a function of these two angles. 

The expressions for the inviscid pressure 
over half-pointed cones are applicable to 
other types of geometry (such as a half 
elliptical cone) if the geometric angles 0 
and 0 at any locations are defined similarly 
to those of the half cone. This implies the 

use of the same axis of symmetry with respect to the centerline. The body local angle 
of attack, \|r, for a semi-cylinder is obtained by setting 9 = 0 in the equation of i|r for a 
semi-cone. 

The thermodynamic properties of the air ahead of the shock, P^, T^, and and the 
free stream Mach number M^, are determined from the given flight conditions using 
the 1963 Patrick AFB atmosphere. 



Figure 2-34. Body Orientation 


a. Low Angle of Attack. A set of empirical equations (Reference 17) is used in 

determining the shock layer thermodynamic properties at the boundary layer edge 
as functions of the hypersonic similarity parameters sin a (Reference 18). 

For the flat plate configuration, a = 0, it is assumed that the free stream conditions 
exist at the boundary layer edge. 


b. High Angle of Attack. At high angle of attack the pressure is calculated by the 
modified Newtonian flow theory, which is expressed as 


where 

Cp 


max 


— — = 1 + Mqq Cp sin 2 t 

P 2 00 ^max 

oo 


.til fi i 1 

y +1 L 0' + 3 ) mJ,J 


(2-24) 


When the attached shock properties break down at high angles of attack, an 
empirical equation based on data of Bertram and Henderson (Ref. 19) is used 
to compute the shock angle as shown in Figure 2-35. 
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Figure 2-35. Effect of Mach Number on Shock Angle for Flat Delta Wings 
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Once the shock angle is known, the conditions downstream of the shock wave can be 
computed from the oblique shock relations using the upstream conditions. The 
governing equations are simply the normal momentum and energy equations 


’2' p l + f>i v In< 1 - x > 


, 2' l l + T V L (1 ' X2) 


(2-25) 

(2-26) 


where x = Pj/p 2 and die continuity equation has been eliminated from the original 
set of equations. The subscripts 1 and 2 denote conditions up and downstream of 
the shock. 


At low freestreem stagnation enthalpies (i < 705 Btu/lb m ), from the perfect gas 
law, the stagnation pressure after an oblique shock, P s , is 


P s_ 

6M 2 sin 2 6 
00 

3.5 

6 
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,2 2 


2 2 

0 

M sin 6+5 


7M sin 6 - 1 


L 00 J 


L 00 J 


At high freestream enthalpies, real gas effects are important, and Mollier chart 
of equilibrium air (Ref. 20) is used. Equations 2-25 and 2-26 are solved iteratively 
together with the equation of state p = p (i, P). 

Aerodynamic Heating 

Aerodynamic heating rates are computed in two combinations: 

Eckert laminar with Eckert turbulent 
Eckert laminar with Spalding-Chi turbulent 


The Eckert reference enthalpy method (Reference 21) depends upon the assumption that 
the incompressible mass, momentum, and energy equations can be used for compressible 
flow solutions provided the thermodynamic and transport properties of the gas are known 
and are evaluated at an appropriate reference enthalpy. Real gas effects, including 
dissociation, are taken into account in the determination of the properties of the gas 
just outside the boundary layer. The Eckert reference enthalpy is 


i* = i e + 0.5 (i w - i e ) +0.22 (i r - i e ) 


e' 


(2-28) 


where i e is the shock layer (just outside the boundary layer) static enthalpy, i w is the 
enthalpy of air at the wall temperature, and i r is the boundary layer recovery enthalpy, 
given by 


ir = ie + r 2ij (2 " 29) 

where the flow recovery factor, r is 0. 84 for laminar flow and 0. 89 for turbulent flow. 
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The shearing stress at the wall is given by 


Cf 2 

r w = - p V e (2-30) 

where Cf is the skin friction coefficient. In the following sections, methods of computing 
Cf will be described. The convective heat transfer is given by 

q = h (i e - i w ) (2-31) 

The heat transfer coefficient is obtained from the Stanton number 


h = p V e St (2-32) 

The Reynolds analogy relation is used to obtain Stanton number 
/ *2/3 

St = Cf/2Pr (2-33) 

where the Prandtl number 
pi* C * 

Pr* = (2-34) 

A constant value of 0. 71 is assumed for Pr*, based on the data of Hansen (Reference 22). 

a. Laminar Boundary Layer . The Blasius solution for an incompressible laminar 
boundary layer gives a skin friction coefficient 

C f = 0.664/ /Re* - (2-35) 


where Reynolds number is 
p* V X 

Re * = —£T~ ( 2 - 36 ) 

yields the laminar heat transfer coefficient 
13.43 yp* pi* V 

h — (2-37) 

JT 


b. Turbulent Boundary Layer. The Schultz -Grunow solution for the local skin friction 
coefficient is (Reference 23) 


7.485 p* V 


C f = 


e 


(log Re* ) 


2.584 


(2-38) 


The Spalding- Chi (Ref. 24) method was based upon incompressible skin friction 
relationships. Transformation from incompressible to compressible cases was 
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postulated upon previously successful transformations of Van Driest and others. 
These transformations were evaluated by obtaining constants with the least squares 
deviation from a wide variety of experimental data. Equations used in this analysis 
stem from curve fits to the Spalding-Chi values developed by Komar (Reference 25). 
In this manner, successive solutions of the equations 




- ( 2 c + b) 

- Jh 2 - 4 ac - 



(2-39) 


where 
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Cf F = exp 


t *) Ai ' l l 
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a *w/*e 


r 2 A w 

b = l+-(y-l) M - — 
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(2-40) 

(2-41) 

(2-42) 

(2-43) 

c = -7 (y-1) M 2 (2-44) 

Li 6 


The coefficients gj are those developed by Komar for the curve fits. 


Reynolds analogy is then used to obtain the heating coefficient. A recommended 
equation for the Reynolds analogy factor is that of K arm an as proposed by Bertram 
(Reference 26). 


S =1 + 5 /(F c c f )/2 


Pr - 1 + log 


5 Pr + 1 


(2-45) 


Cross-flow and Transition Boundary Layer . At vehicle angles of attack greater than 
twenty degrees, for body shapes of the type associated with lifting entry, streamline 
divergence, or cross-flow, thins the boundary layer and increases heat transfer. 

Two dimensional flow assumptions, ie, strip theory, do not account for this three- 
dimensionality and therefore underpredict observed heat transfer. Approximate 
methods have been developed for estimating the influence of flow divergence. These 
methods are described by Young in Reference 27. Since cross-flow also causes 
a reduction in boundary layer momentum thickness, it must be included in any 
transition correlation involving this parameter. 


For a slab of width D, the velocity gradient along the center line is defined by the 
relation 
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The ratios of heating rate with cross-flow to that from strip theory are 
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where 


A =f 0 — 


X 
2 D 


(2-46) 


(2-47) 


(2-48) 


and the ratio of momentum thickness for a slab including cross-flow to that for 
strip theory is 


e CF 

e. 


l-e 


•2A n 1/2 


(2-49) 


An interim transition criterion has been proposed by Masek (Reference 28). The 
value of the transition criterion at which transition onset occurs is a function of flow 
deflection angle as shown in Figure 2-36. The criterion can be expressed as 


Re / M e 

(Re x /X)- 2 


= f (♦) 


(2-50) 


All quantities are based upon local boundary layer edge properties where local 
Reynolds number and momentum thickness Reynolds number are 


and 


Re x 


P V X 
e e 


^e 


Re* 

Re = — 9 

e x 


(2-51) 


(2-52) 


The momentum thickness is given by 

. 664X Ip *11 * e CF 
" (Re ) 1/2 V P e "e ' « s 

X 

where p*, fj,* are evaluated at Eckert’s reference enthalpy. 


(2-53) 
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Figure 2-36. Correlation Parameter 

The transitional boundary layer is represented as a linear transition from laminar 
into turbulent flow. Transition is assumed to begin at a calculated value of the shock 
layer Reynolds number, Ret r , and end at a calculated value of the shock layer Reynolds 
Humber, denoted Reg. The heating rate is computed as the linear interpolation 
between laminar and turbulent values, the virtual origin of both being the leading 
edge. Interpolation is performed using the value of Reynolds number that occurs 
between the calculated values specifying the onset and end of transition. Hence, 
for example, the heat transfer coefficient is given by 

h = Re ~ Re tr + I" Re - Re tr 
Re E ” R e tr ^ L Re E “ Re tr 


h L - f T Itp + fn (2-54) 


The value of Re tr can be determined by Equation 2-50 and experimentally, one 
observes (Reference 29) 

Re E 

f tr = rI“ = 2 (2-55) 

Ke tr 

d. Sphere and Cylinder Heating . For leading edges or at high angle of attack, swept 
cylinder theories are used to predict the aerodynamic heating. For laminar flow, 
heating is computed by method of Fay and Riddel (Reference 30), while the method 
of Beckwith and Gallagher (Reference 31) is used to compute the turbulent heating. 
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Bushnell (Reference 32) suggests the use of a simple transition criterion for leading 
edges at all sweeps which is based upon freestream Reynolds number and the 
leading edge diameter as the characteristic length. Based upon Bushnell's and 
other's data and since fin and wing leading edges of the shuttle will have end effects, 
the onset of transition is taken as R^, D = 0. 1 x 10 6 , and fully turbulent flow occurs 
as Re^, D = 0.2 xlO 6 . 


Laminar heat transfer to a swept cylinder is computed through a transformation of 
spherical heating rates. Heat transfer to a sphere is calculated by the Fay-Riddell 
(Reference 30) expression 


q sph =30.25 (P 8 M S )°' 4 WvW!' 5 U 
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s 


w' 


(2-56) 


where 
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(2-57) 


Cylindrical heating rates are then obtained by adjusting for sweep by 

q = 0. 75 q L (cosA) 1 ’ 2 (2-58) 

cyl sph 

where the sweep angle A = 90° - i|t. 

Turbulent heat transfer to a cylinder is computed by the method of Beckwith and 
Gallagher (Reference 31). The equation for heat transfer to the stagnation line of a 
swept cylinder is 


h D 
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Pr 


1/3 


rp V D-i n/(l+n) r u p -i n/(l+n) t- -i(n-l)/(n+l) 


O 00 


L ^o 
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->n/(l+n) 


(2-59) 


where Pr = 0. 71, a and n are constant in the Blasius skin friction law and were 
taken as a = 0.0228 and n = 4. V n is the component of normal to the cylinder 
stagnation line; thus V n = (cos A ). 


The heat transfer to the stagnation line is then 



= 1.30622 


r V Sin An 3/5 

-2 1 

L ^o J 




(2-60) 


Reference conditions p r and fip are evaluated at stagnation conditions after the shock, 
(Ps, Ts), where values of recovery factor used in evaluating ij. are obtained by curve 
fit of the data in Beckwith and Gallagher (Reference 31). 
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The methods of Bertram and Henderson (Reference 19) are used to correct the flat 
face velocity gradients (Figure 2-37). Since stagnation heating is proportional to 

free stream velocity gradient, it was 
postulated that 



q flat _ 

/ dv e 

/ dx x=0 

flat 

q cyl ^ 

/ dV e 

, e x=0 
dx 

cyl 
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Figure 2-37. Geometry of Flat Face 
Velocity Gradient Correction 


flat 


/ 0.745 + 3.' 
%yl y 2.3 


For the laminar boundary layer of a flat 
disc 
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the laminar centerline heating rates are 
of the form 


rdV -1 
e 

= 0.745 + 3.14 

r 

dx 

x=0 

d” 


14 (r/D) 
315 


(2-63) 


For turbulent boundary layer a 1/5 power correction is used, i.e. , 

1/5 


q flat q cyl 


0.745 + 3.14 (r/D) 
2.315 


(2-64) 


e. Heat Distribution on Peripheral Locations. Lees (Reference 33) solution modified 
by Bertram and Everhart (Reference 34) is used to calculate the distribution of 
heating rates for various peripheral locations. If the peripheral pressure 
distribution over a body is assumed to be represented by an equation of the form 

p/p Q = cos 2 0 + b sin 2 0 + c 0 n (2-65) 


the ratio of heat transfer coefficients at the peripheral location to the stagnation 
line is 


h 

h 

o 


J~2 [1 - (l-b) sin 2 0 + C0 a ] 0 


(l-b)(/3 sin 20 - sin 2 0) + 0 2 ( ~ o 0 n + 1 + b) 


1/2 


(2-66) 


According to the Newtonian approximation adopted by Lees, b = (y M^) and c = 0. 
However, the values of b, c, and n which accurately represented the pressure 
distribution on a cylinder for Mqq s 4 were found to be b = 0. 16, c = -0.00665, 
n = 5 for 0 < 9 = 90°. ... 
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Thermodynamic and Transport Properties 


The perfect gas law is used for air at low temperatures, and the enthalpy of air is 
obtained from an empirical equation as shown in Figure 2-38. At high temperature a 
Mollier chart (Reference 20) is used. 


A constant value of 0. 71 is assumed for Prandtl number, based on the data of Hansen 
(Reference 22). The viscosity of air is calculated by equation of Hansen and Heims 
(Reference 35) which is in the form 




1 + 0.023A 


1+ tanh 


A(l--B)-6.5 

_ 

1.5+-B 

8 


1 - exp 


A-14.5-1.5B 
0.9+ 0. IB 


(2-67) 


where 


A = T/1800 


B = log (P/Po) 


T is temperature in degrees Rankine 


P,Po are local pressure and one atmosphere pressure respectively 


takes Sutherland's value, i.e.. 
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( 2 - 68 ) 


m 3/2 

U = 2. 27 X 1CT 8 — 

Mo u T + 198.6 

2.2 ADDITIONAL TPS DESIGNS 

The capability of treating a number of alternate TPS concepts was added to the existing 
TPS sizing routine. These included the passive systems of Reusable Surface Insulation 
(RSI), the carbon-carbon systems including a new leading edge concept as well as the 
configurations available for the flat panel metallic systems, and a heat sink system 
which incorporates either a change of phase or change of state in a material embedded 
in the honeycomb substrate which supports the RSI. The only active cooling concept 
considered was the charring ablator. 

Of the new systems investigated, the most difficult task by far to accomplish was the 
stress analysis required for the RSI. Because of the large deformation of the glue line, 
a finite element stress analysis had to be shortened and streamlined to allow its use in 
an already massive computer code. In addition, a discrete element stress analysis on 
the same system had to be developed and utilized just to find the design point at which 
the finite element analysis could be used. The next most difficult task was the tempera- 
ture description of the ablator. Fortunately, a charring ablator code was available to 
be adapted into the TPS sizing analysis. Simply by suppressing the chemical reactions 
of the ablator, a one-dimensional heat transfer analysis could be employed for the RSI 
system both with and without the phase change material (PCM). The PCM itself was 
described simply by performing an energy balance on the substrate to determine when 
enough heat had entered this region of the system to support either a change of phase 
or a change of state. Material properties were then adjusted accordingly. Finally, a 
thermodynamic model was developed to describe the carbon-carbon leading edge. The 
following paragraphs describe the mathematical models of the new TPS's. 

2.2.1 PASSIVE COOLING SYSTEMS 

2. 2. 1.1 Reusable Surface Insulation (RSI) 

Thermal Analysis 

No new programming was required to perform the thermal analysis of RSI. This is 
accomplished using the ablator routine by inputting the appropriate RSI thermal 
properties (ie, thermal conductivity, specific heat, etc), and fictitious ablation 
parameters to suppress the chemical reactions and pyrolosis. 

Stress Analysis 

The stress analysis of the RSI TPS is performed by two separate routines. The first, 
already developed within the framework of analyses existing in the TPS sizing driver 
program, describes by discrete element analysis the stress produced within an isolated 
section of a simply supported beam. The loadings considered are bending due to 
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external aerodynamic pressure and internal forces due to temperature gradients. The 
discrete element analysis is performed throughout the trajectory at time increments 
dictated by the program user. A number of structural indices are calculated at each 
time step, and the critical design point of the system is identified as that trajectory 
point exhibiting the minimum design indices. At this time, the second stress analysis 
is performed; here, a finite element analysis adopted from an existing independent 
program is applied. The RSI panel cross-section is subdivided into a group of bars 
and panels simulating the behavior of the RSI, its coating and underlying bond and the 
honeycomb substrate. The advantages of the finite element analysis are that it considers 
axial, bending, shear, and normal deformation of the panel and RSI as well as the shear 
and normal deformation of the bond. However, within the constraints of a fast and 
inexpensive computer program, it is not economically feasible to perform the finite 
element stress analysis throughout the trajectory. Hence, at least initially, this 
procedure will be employed only at the design point of the trajectory. 

During the stress analysis of the program, the following functions are performed. At 
the first call for the stress analysis, subroutine PANEL divides the RSl/bond/substrate 
configuration into a predetermined number of discrete elements. The areas and 
location of these elements are computed as well as the appropriate moments of inertia. 
Once these geometric parameters are determined, there is no need for them to be 
repeated each time step; the values are stored for future use. At each time step, nodal 
temperatures for the discrete element stress analysis are interpolated from the thermo- 
dynamic analysis and stored. Next, thermal and mechanical stresses are computed. 

The stress model is comprised of five different components: (1) the face sheets of the 
honeycomb substrate, (2) the core of the substrate, (3) the bond material, (4) the RSI, 
and (5) the coating which is considered to cover only the top of the TPS. Inherent in 
this assumption is the neglect of edge effects of the coating around the edges (perpendicular 
faces) of the RSI tiles. Material properties are evaluated at the nodal temperatures and 
include the ultimate tension, compression, and shear strengths, the modulus of elasticity, 
the coefficient of thermal expansion, the modulus of rigidity, and Poisson's ratio for the 
bond material. Some modifications had to be made to the existing stress analysis since 
the original dealt with only a single material. A sketch of the nodal breakdown for the 
five different materials is given in Figure 2-39. Once the stresses are computed for the 
nodes of the model, a comparison is made to the applicable allowable values. The design 
factors or margins of safety are presented in Table 2-24. These are computed throughout 
the trajectory and stored for investigation at the trajectory end. The design point is then 
identified by the minimum design factors. 

At this point in the computer program, the finite element stress analysis is performed. 

It is used primarily as a check of the design point which has been located by the more 
economical discrete element method. 

Finite Element Stress Analysis General Theory 

This section presents the general theory for determining displacements and stresses in 
the two-dimensional structure used to simulate the RSI behavior by the finite element 
stress analysis. 
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Figure 2-39. Discrete Element 
Stress Analysis Nodal 
Schematic 


The method is based upon Castigliano's Theorem 
of Equilibrium (Reference 36). Finite differences 
are developed for the partial derivatives, and the 
resulting equations which relate forces and dis- 
placements are written in matrix form for automatic 
digital solution. The program output is internal 
loads and point displacements. 

The development is presented in two sections: 

a. An orthogonal stiffener grid 

b. Thermal stress considerations 

This procedure is advantageous because it 
yields separate routines for the simpler (and 
more commonly encountered) cases, and thus 
saves computing machine time. 

The general development requires four steps: 

1 . Obtain a statement for strain energy at a 
typical point in the structure. 

2. Apply the Energy Theorem. 

3. Express the results of the preceding steps 
in matrix form. 

4. Present the problem in a proper form for 
use by a digital computer. 


The last two steps offer considerable latitude in procedure and have been tailored to 
suit existing computing equipment and programs. 

Development of the Theorem 


Castigliano's Theorem of Equilibrium states, "If the strain energy of any elastic 
structure is expressed in terms of the displacements of a number of points of the 
structure, then the differential coefficient of the expression with respect to the 
displacement of one of the points gives the component force acting at that point in 
the direction of the displacement. " If u and w are displacements in the x and z 
directions at a point (m, n), then 
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Table 2-24. Design Factors 
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It is essential to note that, when using difference equations, the expressions stated 
above indicate that only the structural elements containing or including point (m, n) 
need be considered. (The derivatives of the terms for all elements not containing 
(m, n) are zeros). This permits use of one set of general equations to represent any 
point in the structure. 

Orthogonal Stiffener Grid 

Figure 2-40 shows all structural elements which must be taken into account when a 
typical point (m,n) is considered. 

The following limitations and assumptions apply: 

1. All applied forces are in the plane of the structure. 

2. Stiffeners carry axial load only. 
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Figure 2-40. Structural Elements 

3. The axial load carrying capacity of the panels must be predetermined and included 
as part of the axial load carrying capacity of the stiffeners; the panels carry only 
pure shear. 

4. The stiffness of each axial or shear element is constant for the length of the 
element. 

5. The axial load in a stiffener varies uniformly so that the shear in each attached 
panel is constant. 

6. No stability failure can occur. 

7. There is a constant linear load-deflection relationship for each structural element 
at all times. 

8. Each stiffener is either parallel to or inclined at 90 degrees to each other stiffener. 

For plane stresses in two direction 

_ au # _ 5w - _ B_u + 
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The strain energy expressions needed are 
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For a typical stiffener in the x-direction, for example, bar (m,n) to (m+l,n) 
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and 


U = - C (u - u )* 
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For a typical stiffener in the z-direction, for example, bar (m,n) to (m, n+1). 
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For the shear panel (m, n) 
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When constraints exist at (m, n) the strain energy in the constraints 
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Then for any force applied at point (m, n), the total strain energy involved can be 
expressed as the sum of the energies of the elements surrounding the point, as stated 
previously in the development of the Energy Theorem. 
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All of the terms needed to satisfy these conditions are shown in Figure 2-40. 
It follows from Equations 2-70 and 2-71 
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Equations 2-76 and 2-77 are typical equations for any point in the structure. If the 
point under consideration is on an edge or borders a cut-out, the stiffness parameters, 
C m , n an d n become zero. 

Thermal Effects 


The effects of temperature differentials in a structure can be taken into account for 
an orthogonal grid by a modification of Equations 2-76 and 2-77. It should be noted 
that the thermal effects do not change the stiffness matrix, and hence are in effect 
external forces. 



Ax=Au-9kAx 
Az = Aw - 9 k Az 


9 is the temperature differential of the bar ends, and k is the coefficient of thermal 
expansion. 

U _ = 7 - C u -u - 9k (x -x )"1 

m,n,x,T 2 m,n,x m+l,n m, n m+l,n m,n| 
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Again assuming that the panels carry only shear, their energy expressions do not 
change. They the total strain energy to describe the state at point (m,n) is 
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Then, from Equations 2-76 and 2-77 
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F = F + C T - C ■ T , (2-79) 

m,n,z,T m,n,z m,n,z m,n,z m,n-l,z m,n-l,z 

Formation of the Matrices 


If the unknowns are listed in a specific order, the equations may be expressed in 
matrices of the form 
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Since displacements are the unknown quantities, the problem is solved 



Examination of Equations 2-76 and 2-77 shows that for Point (m, n) a coefficient for 
u and w at each of the nine points considered is found in the expression for F x and for F z . 

After the equations have been solved, the deflections of each point for which Strain- 
Energy equations were written are known. These displacements and the basic 
relationships are used to find internal stiffener loads and panel shear flows. These 
will be shown for the orthogonal grid. 
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This is the internal load at the midpoint of stiffeners (m,n), and a positive sign 
indicates tension in the stiffener. 


For a shear panel 
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These are active shear flows, and the positive sign indicates shear flows in the 
positive direction act on sides with the more positive coordinates. 

Equations 2-79 and 2-80 are valid. 

For the thermal effects considered 
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P = C (w - w - T ) (2-82) 

m,n,z,T m,n,z m,n+l m, n m,n,z 

A special purpose version of the previously described finite element computer program 
(Program No. P2354) has been prepared to perform the stress analysis of a TPS panel 
with bonded RSI. The major modifications to the program consist of (1) the addition of 
subroutines to generate the mathematical model, applied loads and temperatures 
internally from a minimum of input parameters, and (2) the addition of routines to 
compute maximum principal stresses and maximum shear stresses in the RSI and 
bond and to select the critical stresses in each element. Details of the necessary 
input data, mathematical model, stress calculations and an outline of the method of 
incorporating the RSI analysis in the existing stress program are discussed below. 

A comparison of preliminary results with a published analysis is also given. 

The input parameters required to describe the panel geometry are as shown. 

SL = Semi Span (in) TB = Bond Layer Thickness (in) 

TF1 = Inner Face Sheet Thickness (in) TI = Insulation (RSI) Thickness (in) 

TF2 = Outer Face Sheet Thickness (in) TC = Protective Coating Thickness (in) 

H = Honeycomb Core Dept (in) 
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Figure 2-41. Panel Geometry 


The input parameters defining the applied loads are 


PX = Differential Pressure on Panel (psi) 

PZ = Axial Load in Honeycomb Panel at Edge (lb) 




(1) This temperature used 
for "X" bar. 

(2) This temperature used 
for "Z" bar. 


* Z 


Figure 2-42 . Transverse Temperature Distribution 

The input parameters required to describe the transverse temperature distribution 
are 


TX(I) = Temperature at Mid- Length of Each "X" Direction Bar Segment (°F) 
TZ(I) = Temperature at Each "Z" Direction Bar (°F) 

The temperature distribution is assumed constant in the "Z" direction. 
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The input par am eters required to describe the temperature dependent material 
properties are 

NT = Number of reference temperatures at which material properties 
are defined (the maximum number is 16) 

El (I) 1 = 1, NT Elastic Modulus of Face Sheets (psi) 

E2(I) 1 = 1, NT Elastic Modulus of Core (psi) 

E3(I) 1 = 1, NT Elastic Modulus of Bond (psi) 

E4(I) 1 = 1, NT Elastic Modulus of RSI (psi) 

E5(I) 1 = 1, NT Elastic Modulus of Coating (psi) 

G2(I) 1 = 1, NT Shear Modulus of Core (psi) 

G3(I) 1=1, NT Shear Modulus of Bond 

G4(I) 1 = 1, NT Shear Modulus of RSI 

ALFAl(I) I = 1,NT Temperature Coefficient of Face Sheets (in/in-°F) 
ALFA2(I) I = 1,NT Temperature Coefficient of Core (in/in-°F) 

ALFA3(I) 1=1, NT Temperature Coefficient of Bond (in/in-° F) 

ALFA4(I) 1= 1, NT Temperature Coefficient of RSI (in/in-°F) 

ALFA5(I) 1=1, NT Temperature Coefficient of Coating (in/in- ° F) 

VB Poisson's Ratio of Bond 

The mathematical model which represents one-half of the symmetrical panel is shown 
on Figure 2-43. Seven sets of longitudinal bars represent the panel face sheets, RSI 
and the protective coating. The number of sets of transverse bars (NH) is defined at 
input (Min. = 16 Max. = 40). The boundary conditions representing a simply supported 
panel are determined by eight constraints, one at each "Z" direction bar at the panel 
centerline which prevents axial displacement while permitting free displacement in the 
"X" direction. Displacement in the "X" direction is constrained at the panel edge. 

The derivation of the grid spacing, bar areas, bar and panel properties and the applied 
loads from the input parameters follows. 

Distance between sets of "Z" direction bars DX(J) 

DX(1) = H + .5 (TF1 + TF2) 

DX(2) = TB 

DX(3) = DX(4) = DX(5) = 2. * (TI + 0. 5 * TC)/ 7 

DX(6) =0.5*DX(3) 

Distance between sets of "X" direction bars (DZ(I)) 

The basic grid size is progressively reduced as the panel edge is approached. 

N1 = NH - 15 
N2 = N1 + 4 
N3 = N2 + 6 
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DZ(I) »* SL/FLOAT (N1 + 4) 

IF (I.GT.N1) DZ(I) = DZ(l)/2 
IF (I.GT.N2) DZ(I) = DZ(l)/3 
IF (I.GT.N3) DZ(I) = DZ(l)/4 

Areas of "X" direction bars AH(I, J) 

NH1 = NH-1 
AH(1, J) = DZ(l)/2. 

AH(I, J) = .5 * (DZ(I-l) + DZ(2)) I = 2, NH1 J = 1,7 

AH(NH, J) - DZ(NHl)/2. 


Areas of "Z" direction bars AV(I, J) 1 = 1, NH1 

AV(1, 1) = TF1 
AV(1, 2) = TF2 
AV (I, 3) = DX(6) 

AV (1, 4) = AV (I, 5) = AV(I, 6) = DX(3) 

AV(I, 7) = TC 


The temperature dependent material properties required in the station are 


EH(I, J) = Elastic Moduli of "X” direction bars 
CH(I, J) = Temperature coefficient of "X" direction bars 
EV(I, J) = Elastic moduli of "Z" direction bars 
CV(I, J) = Temperature coefficient of "Z" direction bars 
G(I, J) = Moduli of rigidity of shear panels 


Given the temperature distributions TX(I) and TZ(I), the instantaneous values of 
parameters are computed by linear interpolation from the input material properties 
at the reference temperature (T1(I)). 


For the "X" direction bars representing the bond layer (EH(1, 2)), a correction factor 
is applied to account for the fact that under normal loading lateral deformation of the 
relatively soft bond is restrained by the adjacent faces. 


Assuming the lateral deformations are zero the correction factor is given by 


Factor = 


(1 - VB) 

(1 + VB) (1 - 2VB) 


where VB = Poisson's ratio for the bond. 
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Stress Calculations 


The output from program P23 54 consists of the following stresses 

FH(I, J) = Axial stresses in "X" direction bars 
FV(I, J) = Axial stresses in "Z" direction bars 
FS(I, J) = Shear stresses in panels 


For the RSI, maximum principal stresses and maximum shear stresses (Figure 2-44) 
are derived as follows 


1/2 


FP(I, J) = (FV(I, J) + FH(I,J-l)/2. ± 


( 


FV(I, J) - FH(I, J-l) \2 


) 2 + FS(I,J-1) 2 


FSM(I > J)=[( FV ^ J )- 2 FH ^ J - 1 > ) 2 + FS(I,J-1) 2 ] 


1/2 


▲ 


FH (I, J-l) 


FS (I,J-1)- 


► 

V 


-► FV (I, J) 


Due to the repetitive use of the stress 
analysis in the overall TPS optimization 
program it is considered that the sole use 
of the finite element program (P2354) would 
result in prohibitive computing time; hence, 
a basic analysis will be performed using 
the existing thermal stress and overall 
bending analysis routines adapted to the 
RSI panel. Upon completion of the 
trajectory and when the most critical . 
trajectory points (minimum margin of 
safety) are determined, the stress analysis 
at the most critical points will be repeated 
using Program P2354 and the resulting 
stresses used to compute the final margins. A flow chart is given in Figure 2-45. 


▼ 


Figure 2-44. Elemental Volume 


The failure modes considered, the corresponding margins of safety and the buckling 
equations used for the sandwich face sheets are given below. 


Element 


Modes Considered 


Face Sheets 

Core 

Bond 

Insulation 

Coating 


Yielding, ultimate tension, wrinkling, dimpling 
Ultimate shear 

Ultimate shear and normal stress (peeling) 
Ultimate tension, compression, shear 
Ultimate tension, compression 
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Figure 2-45. Flow Chart of RSI Analysis 
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Face Sheets 
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F NET 
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where FCR is given by the smaller of the critical wrinkling or dimpling buckling 
stresses. 
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- 1 
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The sandwich panel face sheets are analyzed for two modes of local buckling under 
compression loads as follows: 

a. Wrinkling. In this mode the face sheets buckle in a series of sinusoidal waves. 
The critical buckling stress is given by 

V " °‘ 5 <°C E cz V E f ) 1/3 

where 


G c = modulus of rigidity of the core 
E cz = modulus of elasticity of the core 
Ef = modulus of elasticity of face sheet 
77 plasticity factor = E-p/E 
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b. Ihtercell Buckling. In this mode the face sheet distorts in a pattern of dimples the 
size of the honeycomb cells . 


The critical buckling stress is given by 
a cr = 3,1 ? E f (t f /d > 2 


where 

d = honeycomb cell diameter 
tf = thickness of face sheet 
Ej and rj as given for wrinkling 

To evaluate the finite element model of Program P2354, a comparison was made with 
a published analysis of Lockheed, (Reference 10). 


Verification of Finite Element Stress Analysis 


Preliminary runs using Program P2354 have been made simulating the model used on 
page 3-3 of the reference as shown below 


80,000 LB /IN 



The structural panel used in the model is represented by a solid 0.125 inch plate. 
Therefore, this plate must be idealized as a sandwich panel for Program P2354. 

For combined bending and axial loading, the bending and axial stiffnesses must be 
represented correctly. In addition, the distance from the bond line to the panel 
centroid must be preserved in the idealized model. 
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The requirements are: 


(1) 

•^1 + ^2 ^P 

(2) 

A 2 (Y) 2 +A 1 (H-Y) 2 =Ip 


A 1 H A 1 H - 

(3) 

- Y 

A. +A„ A_ 


12 P 

This yields 


H„ = 


2 "AjY 


+ Y 


Y • A 

a i = _ h7" 


a 2 = A P (i - Y/H) 
For 0. 125 Plate 



Ap = 0.125 in 2 

H = (°- 1253 / 12 ) 

c 0.125x0.0625 


Y = 0.0625 in 
+ 0.0625 = 0.08333 in 


A 


1 “ 


0.0625 xQ. 125 
0.08333 


= 0.09375 in 


A 2 = 0.125 (1 - 0.0625/0.08333) = 0.03125 in 

The shear stiffness of the idealized honeycomb panel is represented by using an 
equivalent modulus of rigidity for the core . 

G _ ° P Ap - °- 125 G 
e H 0.0833 p 
c 

The overall agreement is good except for the bond normal stress (approximately 60% 
of Reference 10 values). This is considered due to the fact that the P2354 analysis 
did not consider lateral restraint of the soft bond layer (E = 300 psi) by adjacent faces. 

Assuming lateral deformations are restrained, the effective modulus of elasticity is 
given by (Reference 37) 

E - E a - r) 

ff (l+y)(l-2 y) 
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Table 2-25. Comparison of Maximum Bond and RSI Stresses (Neglecting Correction 
to Bond Elastic Modulus Due to Lateral Restraints) 


Case 

No. 

Bond 

Thickness 

Bond Stresses 

RSI 

Stresses 

Max. Shear 

Max. Normal 

Max. Shear W 

Max. Principal 

Al 

0.032 

40.1 (41) 

26.4 (49) 

40.1 (38.0) 

122 (111) 

A2 

0.020 

51.3 (51) 

27.1 (60) 

51.3 (42.0) 

139 (125) 

A4 

0.015 

59.0 (59) 

45.1 (70) 

59.0 (56.0) 

147 (138) 

Values in parentheses are from Reference 3 7. 



(1) The maximum shear in the RSI occurs at the bond line and is equal to the bond 

shear stress. 





Table 2-26. Comparison of Maximum Bond and RSI Stresses (With = E e ££ 

and=0.45) 


Case 

No. 

Bond 

Thickness 

Bond Stresses ^ 

RSI Stresses 

Max. Shear 

Max. Normal 

Max. Shear ^ 

Max. Principal 

Al 

A2 

A4 

With t 
as sun 
effect 

0.032 

0.020 

0.015 

his correct 
ted in order 
on the othe 

41.2 (41) 

52.8 (51) 

60.8 (59) 

ion a value of J 
to match the [ 
r peak stresse 

41.6 (49) 

57.7 (60) 
69.2 (70) 

Poisson's ratio 
janel bond stres 
s is negligible. 

41.2 (38) 

52.8 (42) 

60.8 (56) 

approximately e( 
ses given by Rei 

124 (111) 

140 (125) 

148 (138) 

}ual to 0 . 45 must be 
ference 37. The 


where 


E = bond modulus of elasticity 

y = Poisson's ratio 

2.2. 1.2 Carbon-Carbon Composites . The carbon-carbon leading edge was investigated 
for adaptation to the TPS sizing routine. Carbon and graphite composites are important 
because they possess excellent high temperature stability and retain structural integrity 
at temperatures above 4,000F. This latter property alone makes the carbon-graphite 
class of unique candidated TPS materials. To fully utilize the high temperature 
properties of carbon and graphite composites for reusable entry thermal protection 
systems, an impervious oxidation barrier must be provided for the material. Efforts 
have been made in the past few years to increase the oxidation resistance of carbon and 
graphite composites by the use of select fillers and coatings. These effects have been 
successful in providing improvements in oxidation protection. To provide the capability 
for multiple, long-time exposures to expected shuttle entry environments, advanced 
oxidation protection techniques are currently being developed. 
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Because of the problems encountered with laminated carbon-carbon composite 
material (limitations in design flexibility associated with the complex tooling, local 
reinforcements for attachments, and overall reliability) a 3-D reinforced carbon- ' 
carbon fluted core concept has been developed (Figure 2-46) as a lightweight, double 
faced, truss core structure that integrates an external heat shield cover panel with 
a support structure. The support can in turn be attached to a primary load carrying 
structure. The concept consists of a 3-D woven truss core fabric that is rigidized 
with high-carbon-yield resins, pyrolyzed, and then carbonized up to 4, 000F in an 
inert atmosphere. A simulated shuttle vehicle leading edge demonstration article 
has been fabricated to demonstrate the feasibility of forming the truss core. The 
development of the oxidation inhibited carbon-carbon 3-D fluted core has been 
accomplished in a joint effort between the Convair Aerospace Division of General 
Dynamics and the Defense Products Division of HITCO. 

The truss core configuration shown in Figure 2-46 was selected as the baseline 
concept for the TPS optimization studies for two reasons: the truss configuration is 
structurally very efficient, and experience in fabrication has shown it is easier to 
process than other concepts. The thicknesses and other dimensional variables are 
independent, allowing for considerable flexibility in design. However, it is 
emphasized here that almost any carbon-carbon configuration can be analyzed by the 
TPS sizing routine by using the appropriate equivalent properties; the analysis is not 
restricted to truss core. 

Thermal Analysis 

To treat the heat transfer as a two-dimensional problem, the panel is considered as 
three layers; an outer and an inner layer of carbon-carbon material with an "equivalent” 
material between them. The middle layer has the same mass as the truss core which 
it represents. This defines a density and specific heat of the equivalent material. 

The resistance to heat flow is given by an "apparent" thermal conductivity. The truss 
core apparent conductivity and other properties will be modified by geometry changes, 
and by the use of filler materials (for radiation barriers) in the gaps. 
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The basic problems to be treated by the thermal analysis and insulation sizing are: 

1. Compute temperatures of the panel, underlying insulation and internal 
structure: check that panel temperature is below the temperature constraint. 

2. Increase thickness of insulation between panel and internal structure (if 
necessary) to reduce structural temperature to below the maximum 
allowable value. 

The geometric representation of the outer surface airfoil contour is given in Figure 
2-47, and the nodal representation of the panel for the conduction matrix is given in 
Figure 2-48. 

A thermal analysis has been made for 
an example configuration with geometry 
as given in Figure 2-49. This corres- 
ponds to the leading edge of an NACA 
0012-64 airfoil which is taken as a 
typical wing section for a Space Shuttle 
orbiter. 

Trajectory data were taken for a typical 
30° angle of attack entry. The leading 
edge stagnation line heating was 
computed for a laminar swept cylinder 
using the method of Kemp and Riddell, 
and is shown for this trajectory in 
Figure 2-50. The off- stagnation- line 
heating was obtained from the distri- 
bution of Lees shown in Figure 2-51. 

Figure 2-52 shows the temperature 
versus time of three outer panel 
nodes and the temperature response 
of the structural node achieving the 
highe st te mpe r atur e . The se re suit s 
illustrate the temperature gradients 
in the directions along and normal 
Figure 2-47. Leading Edge Outer Surface to the panel surface through to the 
Contour Geometry internal structure. It is of interest 

to note that the results obtained from 

the TPS sizing routine and shown in Figure 2-52 using a calculation time interval of 
0.5 second compared very closely with computations carried out using the General 
Dynamics Variable Boundary II Program (an implicit method of solution). The explicit 
solution of the sizing routine used about half of the computer time of the other program. 
In addition, close agreement was obtained with repeat runs using time steps up to 1.25 
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PANEL CONFIGURATION CONDUCTION MATRIX 

NSEC = 7 



Figure 2-48. Carbon-Carbon Leading Edge Conduction Matrix 


Carbon-Carbon Leading Edge Example 



A. 

Outer Surface Contour Geometry 





R = 5.30 in 
Lj = 3.06 in 
L 2 = 6. 31 in 

0 0 

0 0 

CD CD 

II II 

«-< <M 
a> cd 



B. 

Conduction Matrix 






Dimension in Inches 


Column No. 1 

_2_ _3_ 

_4_ 

5 6 7 8 


fiX 3.25 

3.06 1.85 

1. 85 

1.85 2.77 2.77 3.06 


AYi (C-C Outer Panel) .10 

.10 .10 

.10 

.10 .10 .10 .10 


A Y 2 (C-C Truss-Core) . 46 

.46 .46 

.46 

.46 .46 .46 .46 


A Y 3 (C-C Inner Panel) . 04 

.04 .04 

.04 

.04 . 04 . 04 . 04 


A Y 4 Pynaflex Insulation) 1.00 

1.00 1.00 

1.00 

1.00 1.00 1.00 1.00 


A Y 5 (Tj Structure) . 25 

.25 .25 

.25 

.25 .25 .25 .25 


Figure 2-49. Dimensions for Carbon-Carbon Leading Edges Example 
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Figure 2-52. Reentry Temperature Response of Leading Edge 



seconds; for an input calculation interval of 2.5 seconds, the solution diverged and the 
program automatically adjusted to 1.25 seconds. 

The temperature response of the interior material normal to the stagnation line node 
(nodes 4, 12, 20, 28, 36) at various times during reentry is given in Figure 2-53. 
Figure 2-54 shows the effect on structure temperature of increasing the insulation 
thickness from 1 inch to 2 inches. 

Stress Analysis 

A stress model was also developed but never implemented to the computer program 
due to increased emphasis directed toward the charring ablator and RSI. It is included 
here for the sake of completeness. 

A section of 60° truss core carbon-carbon leading edge at a given span location is sized 
for combined airloads and thermal loads. A single value of facesheet thickness and a 
single value of core truss thickness are calculated for a constant sandwich depth. 

The shape of a chordwise section 
of leading edge is idealized to a 
constant radius section and two 
straight sections (Figure 2-55). 
Figure 2-56 illustrates the 
procedure. The model is then 
divided into a fixed number of 
elements : 5 elements for the 
curved section and 3 elements 
for each of the straight sections. 
Element length, centroid location, 
and angular orientation are 
calculated based on the input basic 
geometry. Upper and lower edges 
are assumed to be fully fixed, 
and the overall sandwich depth 
remains constant as input. 

Section properties for a unit 
spanwise width of leading edge 
are calculated based on an input 
estimate of facesheet thickness. 


Airload pressures acting on each 
of the elements are calculated as 
part of the aerothermal environment 
for each load condition. The 
resulting normal loads are computed 

Figure 2-53. Temperature Response Normal to and resolved into horizontal and 
Surface Node 4 vertical components. 
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TEMPERATURE 




Figure 2-55. Schematic of Carbon-Carbon Stress Analysis Nodal Geometry 
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Figure 2-56. Stress Routine for Carbon-Carbon Leading Edge 


















Element temperature data, including inner and outer facesheet temperatures, is 
computed throughout the trajectory as part of the thermodynamic analysis. Tempera- 
ture changes are found using an input stress-free temperature. Material properties 
are calculated from a table lookup as a function of temperature, and properties for a 
particular element are selected based on the temperature of the outer facesheet of the 
element. 


The redundant fixed-end reactions due to airloads are found using the principles of 
statics and virtual forces. The reactions (forces: Rx U , R zu , R^, and R zL , and 
moments: M^, My L ) at the upper and lower ends of the leading edge model are 
solved from six simultaneous equations 

i=ll 

\n +R KJ, + 1 P ix = 0 (2-83) 

i=l 

i=ll 

P iz = 0 (2-84) 

i=l 


R zu + R 


zL 


I 


M + M , + 
yu yL 



(2-85) 



M 


yL 


+R _Z.-R X. - y 

xL l zL l i-, 


j=*-ir 


5=1 L 


P. (Z.-Z.)H-P 
jx v 1 J jz 






= 0 
( 2 - 86 ) 


2 


i=ll> 


i=l 


v 5=i-l 

- 

Z.L. 

M +R Z -R x.- > 
yL xL l zL l 

3=1 

P. (Z.-Z.)+P (x,-x ) 

jx v l 3 jz v 1 i 

1 1 
E. 

l 



C( 


M +R Z.-R Y 

yL xL i zL 





Thermal loads are found by calculating the unrestrained thermal deflections arui then 
applying the principle of virtual work. 


6 

xu 


i=U -VjtAT^+AT^) 
' 2 


sin 9. 


(2-89) 
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zu 


i=ll L. a. (AT + AT_.) 

V 1 1 Ai Bi 

= > cos 8. 

L,. „ 2 1 

1=1 


(2-90) 


>u=2 

i=l 


i=ll L. a. (AT - AT ) 
l l ' Ai Br 


(2-91) 


Combining the airloads and thermal loads, the internal load distribution is obtained. 


P. = 
i 


M.= 


V i = 


J=i-1 


< R xL +R xLT )Sine i- (R zL +R zLT )C0 Se i + Z. , < V" ^ V° S ^ 


3=1 


(2-92) 


j=i-l 


M 


■yL +M yLT + <R xL +R xLT> Z i-‘ B *L +R zLT» 5( l + 1 (P Jx Z fVl> 


(2-93) 


j=i-l 


(R _ +R Tr Jcos8.-(R T +R ) sin 0. + S (P cos 0 - P sine ) 
v xL xLT i ' zL zLT 1 jx 1 jz 1 


(2-94) 

The facesheet thickness is sized for combined axial load and bending on the sandwich 
section. The core material is not included in the section property and sizing calculations 
for facesheet thickness. Thicknesses are sized for all elements and load conditions 
considering the following failure modes: tension, block compression and compression 
buckling. The largest value of skin thickness is selected as typical of the inner and 
outer facesheets for all the elements at the span location being considered. 

For each load condition called during the trajectory, calculate 


Pi Mi 

N,. = — + 


li 2 0.884 h 

P,- 


n 2 ,= 


Mi 


i 2 0.884 h 


(strength) 


t 3 

t li F. 


N o. 

t = — 
% F. 


(2-95) 

(2-96) 

(2-97) 


where 
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Ftu 11. ifNis < + ) 

F i = 

l F eu Ui U-NlsH 

F T u lli Fcu ll‘ are ^ ase< ^ on ^Ai* 

(buckling) 

For all values of Nj which are negative, calculate 
p 0.380 h 2 |N | -| 1/3 

t3i L“^ J 


(2-98) 


The thickness is then calculated as the largest of either t^, t 2 i , or t 3i . 

The sheet thickness of the truss core diagonals is calculated considering shear and 
crushing due to bending . Thicknesses are sized for all elements and load conditions, 
and the largest value is selected as typical for all elements at the given span location. 

For each load condition, calculate 



(2-99) 

(2-100) 

(2-101) 


<h - 2t f ) w 0 . 

° 2 i 1. 50 F eun. 

The core thickness is then chosen as t c ^ or t c ^ . , whichever is greater. 


2.2. 1.3 Passive Cooling Systems - Phase Change Material . The primary task 
performed during this investigation on the phase change material (PCM) system 
was the thermodynamic analysis of the phase change. No stress analysis was 
indicated or performed. 
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An analytical model has been developed for predicting the transient one-dimensional 
thermal performance of a thermal protection system utilizing a phase change material 
(PCM) in the backup structure. The heat protection system is considered to consist of 
either an ablating material or a reusable surface insulation (RSI) and a backup 
structure or substrate. The backup material is considered to consist of at least two 
face sheets with a PCM in between. Also, a honeycomb filler material (aluminum or 
titanium) can be used to contain and in some cases to improve the thermal conductivity 
of the PCM. 

The computer subroutine used is a modified version of Curry's Charring Ablator 
Program described in Reference 3. This modified form of the program (P5570) can 
handle one ablating or non- ablating material (RSI) and a maximum of 6 backup materials 
including the bond. Conduction, radiation and/or convection are allowed between the 
materials except for the PCM; only conduction is allowed into and out of the PCM. hi 
addition, the program specifies that there must be at least one material behind and one 
material before the PCM in the backup structure. The thermal properties of all 
materials in the backup structure are temperature dependent. 

The basic governing differential equations for the charring ablator analysis, their 
implicit finite difference formulation, the program input and output, the FORTRAN 
program statements and nomenclature are presented in Reference 3. The present 
discussion contains the analytical expressions required to handle the latent heat of 
fusion or transition (i.e. , for a solid/solid phase change material) of the PCM, 
calculation of the equivalent thermal properties of PCM and filler material and the 
changes required in the mathematical procedure to handle the known phase change 
temperature in the backup structure . 

The set of algebraic equations resulting from the implicit finite difference formulation 
of the governing partial differential equations and boundary conditions is solved by the 
Gauss' elimination method specialized for a tridiagonal matrix. 

Figure 2-57 is a schematic of the thermal protection system considered in this analysis. 
The ablating surface is used in conjunction with a backup structure which consists of 
three backup materials in the present analysis. It should be noted that the entire system 
may be composed of noncharring materials. 

The backup structure may consist of up to six different materials (with or without air 
gaps) with a maximum of nine nodes per material. The PCM constitutes a special 
kind of backup material for which heat transfer in and out can only occur by conduction. 
The PCM itself must contain at least three nodes; if the program user wishes to 
simulate face sheets on the PCM, each of these too must contain a minimum of three 
nodes. 

The heat balance equations for a typical node in the backup structure is given by the 
one-dimensional Fourier conduction (neglecting any heat generation term). Thus 
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CHAR LAYER 


REACTION ZONE 


VIRGIN PLASTIC 


PCM 


INTERIOR 

ENVIRONMENT 


Original Material Surface 


T 


T 

S 



i = 1. j - 1 
i = 2, j = 1 


S 



(In this case: 
NPC = 2) 


i = NP, j = 1 
i = NP + NPM(l), j = 1 

j = NP+NPH(1) 
+NPM(NPC), j+NF 
i = NP + NPBS 
j = NBM 


Figure 2-57. Schematic Diagram of Charring Ablator Thermal Protection 
System With PCM in the Backup Structure 



(2-103) 


The implicit difference form of Equation 2-103 for the i th node is 

T» , - T» T! - T» Ax (T» - T ) 

l-l i i l+l „ ' i i' 

- = p C 

Ax Ax Ax Ax P A0 

2k. , 2k. 2k. + 2k, „ 

i-l i i i+1 

where the prime superscript denotes values at the end of the time step. 


(2-104) 


A0 = 0' -e (2-105) 

Corresponding equations for each node in the system are solved simultaneously to 
calculate the temperature at each node. 

Rearranging Equation 2-104 to put the unknown temperatures on the left hand side 
gives 
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(2-106) 


1 


1 

1 

p i c Pi i 

Ax Ax 

Ax Ax Ax 

Ax 

+ 

AS 

-2k. , + 2k, - 


-2 k. , + 2k. 2k. 

2k. . 

- 


Ax Ax 
+ - 


r P t C p . Ax 


i+1 


L2k. 2k. -I 
l i+1 


A0 


T i 


which is of the form 

at! 1 + B t! + C Tj = D (2-107) 

The computer program generates a tridiagonal matrix for the nodes of the system, e.g., 

»Di 


B 1 T 1 + C 1 T 2 
a 2 t 1 + b 2 t 2 + c 2 t 3 


A 3 T 2 + B 3 T 3 + C 3 T 4 


— Dr 


= Dr 
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A N T N-1 + B N T N " D N 


and solves the simultaneous equations. 

Since the PCM is specified to consist of 3 nodes, the material is lumped into two half- 
nodes on each face and one complete interior node. The thermodynamic properties and 
latent heat of fusion or transition of the PCM materials can be obtained from References 
38 and 13. The equivalent thermodynamic properties of the PCM and filler material is 
computed according to the method given in Reference 38. Thus 


k eq = k PCM + (k F “ k PCM ) F FR < 2 " 109 ) 

Cp eq ~ Cp PCM + < C PF " C PPCM^ F FR ; (2-110) 

P eq “ P PCM + ^ P F ” P PCM* F FR (2-111) 

where Fp^ is the filler fraction given by 

F FR = a f /A T (2-112) 


Inherent in these equations are the following assumptions: (1) the contact resistance 
between filler material and the face plates is negligible , and (2) three-dimensional heat 
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transfer effects are neglected, based on the assumption that the filler materials are so 
closely spaced that these effects are negligible. The errors introduced by these 
assumptions are discussed in Reference 13. 



k i,j+l / AT \ 

i+l.j+l \AX/ 


Figure 2-58. Nodal Schematic of PCM 


When the temperature of each 
node in the PCM reaches the 
melt or transition temperature, 
it is assumed that the latent 
heat of fusion or transition 
(solid/solid phase change material) 
is absorbed by the individual nodes 
at constant temperature equal to 
the melt temperature . 

The physical model for the last 
node of the front face sheet and 
first node in the backup structure 
is shown in Figure 2-58. 


For this interface 

T i,j = T i,j+1 = Tmelt < 2 " 113) 

where TjvielT is the melt temperature of the PCM. 


The heat balance equation for this node during the phase-change process is 


T* . - T*. . t' - T' AQ 

i-l. J i.J i.j+l i+l, j+1 pc 


Ac. 

1 


Ax 

+ — L_ 

2k. ' . 2k. . 

1-1,1 i,l 


Ax 


Ax 


A0 


(2-114) 


2k. . „ 2k. . , 

1,1+1 i,i+l 


where i is the nodal index of the material. 


Rearranging Equation 2-114 and noting Equation 2-113 gives 
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where AQp C is the heat fusion absorbed by the PGM material in the node during the 
time interval A9. Equation 2-3.15 also applies for the last node of the PCM and the 
first node of the back face sheet. For the interior node of the PCM a similar analysis 
leads to 
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Ax. Ax. 
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Since the melt temperature is known, one column of coefficients A, B, C in the tri- 
diagonal matrix becomes zero when the PCM starts to melt. Typically, if T 4 «TMELT» 
we have 

B 1 T x + C 1 T 2 =D 1 

A 2 T 1 + B 2 T 2 + C 2 T 3 = d 2 

A 3 T 2 +B 3 T 3 +0 = D 3 

a 4 t 3 + 0 + C 4 t 5 = D 4 (2-117) 

0+B 5 T 5 + O 5 T 6 =D 5 
• • • 

▼ : 

a n t n-i + b n t n =d n 

Since one column is all zeros, the overall matrix contains a singularity and this matrix 
cannot be solved the usual way. To overcome the problem, the matrix can be split into 
two matrices and the Gauss Elimination Method can then be applied to each submatrix. 

For the particular example shown in Equation 2-117, Tg would be obtained from 
solution of the upper submatrix, Tg from the lower submatrix and T 4 is already known. 

2.2.2 ACTIVE COOLING SYSTEM - ABLATOR 

2.2.2. 1 Stress Analysis . No new stress analysis was performed for the ablator system. 
The substrate is sized to provide the total strength requirement of the ablator/substrate 
combination to satisfy the condition at the end of entry where the ablator has been 
theoretically consumed. The stress analysis already existing in the computer program 
is used for this. 
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2 . 2 . 2 . 2 Thermal Analysis . The thermodynamic analysis of the ablating TPS is 
performed by the one-dimensional, implicit computer code developed by Curry and 
described in detail in the technical note NASA TN D-3150 (Reference 3). The finite 
difference approximation to the energy equation is written in backward difference form 
to allow simultaneous solution of the temperature distribution at the new time step in 
terms of those at the previous time. The purpose of employing this technique is to 
overcome the difficulties of numerical instability since the method is stable for all 
time steps. In the original NASA computer code, heating rates were input using cold 
wall values uncorrected for heat blockage due to the addition of pyrolysis gases to the 
boundary layer. This cold wall rate is expressed as 

q=Mf 0 -i3oo) (2-H8) 

where the enthalpy i Q is a fictional total enthalpy, the sum of the enthalpy of air at 300K 
and the freestream kinetic energy. The hot wall convective heating rate is then determined 
by multiplying the above expression by the fraction 


l - l 

o w 

l o - hoo 


(2-119) 


where the term i w is just the instantaneous wall enthalpy. Such a procedure assumes that 
the convective heating rate coefficient is only a weak function of temperature. Two 
additional options have been added during the course of these investigations. The first 
allows the input of hot wall heating rates thus bypassing the correction from cold wall 
values and the need to input vehicle velocity as a function of time. The second option is 
simply the use of hot wall aerodynamic heating rates predicted by the program itself. 


In any case, the hot wall heating rates finally determined by the program must also be 
corrected for the blockage effect of both the exhausting pyrolysis gases and charring 
ablator as well as for radiation to deep space. The latter correction is simply effected 
by subtracting off the term a c T^ w . However, the blowing effect correction incorporated 
into the TPS sizing routine is different from that originally developed by Curry, anH for 
this reason it will be described here in some detail. 


The method for predicting blowing effectiveness published in NASA TN D-3150 is 

Vock^g < r o- (2-120) 

where the driving potential is the difference between the fictional total enthalpy and the 
wall enthalpy. The heat blockage is also a function of the mass flow of the released 
pyrolysis gases and a coefficient tj which has a wide range of possible values. Since the 
heat blockage term is subtracted from the incident heating rate, a high enough value of 
V can e fl®ct a negative net heating rate. Unfortunately, little information is available 
concerning realistic values of the blowing effectiveness coefficient as used in the 
context of Equation 2-120. 


\ 
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Therefore, another blowing effective relationship has been incorporated after a suggestion 
of Curry. The expression is 


S T | 

— = ( 2 - 121 ) 

To e *-l 


where Srp is the Stanton number, the subscript ’’o" denotes the no-blowing case, and the 
parameter § is given by 


5= f 


BLOW 


28.97 j8 
M - 


w 



( 2 - 122 ) 


The parameters fBLOW fi are given in the following table for laminar or turbulent 
flow. 


Laminar Turbulent 

f BLOW 1,44 °' 87 

0 0.25 0.20 


The term M is the molecular weight of the mixture 
of pyrolysis gases and char release to the boundary 
layer, and the product p w v w denotes the mass flow 
released. 


An immediate advantage observed in the relationship of Equation 2-.121 is that the ratio 
of Stanton numbers for the blowing and no blowing case now varies between zero and 
one, the former limit corresponding to massive blowing and the latter to the no blowing 
case. 


In order to include the rigorous approach of the charring ablator code in describing 
ablator performance yet keep the computer program core size and operational speed 
within realistic limits, it was necessary to decrease the size of the ablation subroutines 
and eliminate a number of computations . Hence, the original scheme which allowed up 
to 50 nodes in the ablation material and 12 backup structure materials with 10 nodes 
each was reduced to one which limited the ablator to 10 nodes with six backup materials 
of up to 4 nodes each. In addition, a number of other options of the original program 
were also eliminated. These include the plot routines, the location of constant 
temperature lines as a function of time, the determination of temperatures at discrete 
points in the model to correlate ablator behavior with measured thermocouple data, and 
the input of a linear initial temperature distribution. This last capability is now li mi ted 
to either a constant temperature distribution or specific input values for each node of the 
model. 


Taken altogether, the computer code simplifications and reduction of array sizes 
resulted in a reduction in core size from 65.6Kto 32. 5K octal. Program accuracy or 
speed was unaffected by these changes. A number of cases were run over a 300 second 
trajectory to assess the effects of these changes. The resulting temperatures of the 
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surface node and the backup structures are shown in Figure 2-59. Little difference can 
be seen in the surface temperature for the three nodal spacings examined and the 
temperatures predicted for the backup structure are identical. The time at which the 
surface node disappears is the same for all cases investigated. The depth of surface 
recession is given in Figure 2-60. Again, little difference is seen between models. 
However, in Figure 2-61, the gas ablation rate (the rate at which pyrolysis gases are 
generated) varies markedly between the models when using the original blowing 
effectiveness model presented in NASA TN D-3150. The time increment of the 
computation was also varied in an attempt to stabilize the solution; the predicted gas 
ablation rate is not converging as either the nodal geometry or the time increment is 
decreased. The new blowing effectiveness equations were input to the main subroutine 
of the charring ablator program. The cases of the convergence study were rerun using 
the new method, but oscillations in the generation rate of pyrolysis gases and the net 
heating were still observed (Figures 2-62 and 2-63). Nevertheless, the temperature 
response of the ablator for the variation in the number of nodes (Figures 2-64 and 2-65) 
do not vary from each other catastrophically, especially in the backup material. Hence, 
the model of reduced size (e.g., a maximum of ten nodes in the ablator material) was 
adapted to P5490B. 
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Figure 2-59. Convergence Study of Ablator Temperatures 
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Figure 2-60. Convergence Study of Surface Recession Depth 
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Figure 2-61. Convergence Study of Gas Ablation Rate 
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Figure 2-62. Convergence Study of Ablator Net Heating 
Rate for Laminar Blowing 
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Figure 2-63. Convergence Study of Ablator Net Heating Rate 
for Turbulent Blowing 
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Figure 2-64. Convergence Study of Ablator Temperatures for Laminar Blowing 
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Figure 2-65. Convergence Study of Ablator Temperatures for Turbulent Blowing 
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SECTION 3 


PROGRAM APPLICATIONS 


3.1 SHORT TERM IMPROVEMENTS 

During preliminary sensitivity and optimization runs on the metallic TPS cover panels, 
consideration was given to developing automatic methods to speed the overall assimi- 
lation of TPS design data and to ease the program user's task in generating this 
information. Sample cases were run varying panel and insulation thicknesses and 
panel sizes. The results of these studies included system performance parameters 
such as stress and temperatures, panel weight, manufacturing cost, and total program 
development cost. For a given panel material, configuration, and location on the 
vehicle, the principal contributor to system cost and weight is its size. The optimum 
system was chosen as the one with either minimum weight or cost. Usually, with a 
little foresight and/or experience on the program user's part, this minimum could be 
established with three sample cases, (i.e., three test points). These points were then 
plotted or sketched quickly by hand, and if the minimum could not be seen, the 
concavity of the plotted results indicated at which end of the data (i.e., smaller or 
larger panel sizes) additional designs had to be run. Similar results were also 
experienced when sizing insulation or panel thicknesses to meet temperature or stress 

constraints. Experience at running the 
program indicated that it was far easier to 
determine the minimum insulation thickness 
for a prescribed temperature constraint 
simply by plotting the temperature of interest 
as a function of material thickness (Figure 3-1). 

Techniques which were considered to replace 
the manual decision-making process included 
automatic plotting of the results, curve fitting 
results and differentiating the resultant curve 
to ascertain the minimum, or adaptation of an 
iterative procedure to determine the minimum 
value of the driving parameter. The first of 
these methods, automatic plotting, was rejected 
because such procedures would limit the sizing 
routines to particular computing systems, but, 
more important, was the delay of up to 24 hours 
caused by waiting for results to be plotted and 
reproduced photostatically. Another alternative 
is also available, on most systems, use of a 



Figure 3-1. Example of Primary 

Structure Temperature Versus 
Insulation Thickness 
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cathode ray tube, but for the machine on which the routine was developed (the CDC 
6400) core requirements of 40K octal for such a capability would require such drastic 
overlaying of the program as to make its use inefficient. In addition, there might be 
a tendency for the user to make rather rapid judgments in an effort to run the next 
point and size the system thus leading to errors in judgment and added machine time 
and cost. It is felt that a little contemplation in deciding upon additional test cases to 
be run should the initial values prove insufficient is a more efficient method than 
running a large number of cases to determine one optimum design. 

The concept of curve fitting results to find the design point was also considered. For 
the case of determining the minimum insulation thickness to satisfy a given temperature 
constraint (Figure 3-1) the three data points must be curve fit, and the intersection of 
this curve and the constraint line found either graphically or numerically . To find an 
optimum panel size, Figure 3-2, the curve which fits the data must be differentiated 
to find its minimum. In either case, extrapolation of the results is dangerous at best 


o 
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if the solution is not bracketed. . If convergence of the procedure is not successful, 
a great deal of machine time can be wasted getting additional data points. 

An iterative procedure offers all the disadvantages of curve fitting or graphically 
representing data. The expense of running enough cases to establish the direction of 
the iterative procedure precludes economical use of the program. 

Thus it was concluded that the most efficient manner to size the TPS was to take 
advantage of the man-machine interplay by allowing the program user to exercise 
his discretionary judgment in running point design. 

3.2 FORMAL OPTIMIZATION PROCEDURES 

As a result of studying methods to automatically size the TPS, it was concluded that 
a large number of cases would have to be run to allow the proper system evaluation 
by any formal optimization procedure. Since such methods prevented economical 
program use, this portion of the study effort was eliminated by mutual consent of the 
contractor and monitor. 

3.3 OPTIMIZATION STUDIES 

In the previous section, the investigation of formalized optimization procedures was 
discussed. It was concluded at the end of this study that the wide variety of parameters 
affecting the performance of the TPS precluded the use of anything but the simplest of 
procedures. The basic technique recommended at conclusion of the last contract study 
was to vary the panel thickness to adequately withstand both the mechanical and thermal 
stresses and to vary the insulation thickness to satisfy temperature constraints. All 
other variables such as panel size, configuration, materials, and the like were 
maintained at constant value. The sizing of the panel and insulation thicknesses to 
satisfy temperature and structural constraints was done automatically simply by 
adding to material thicknesses until the constraints were met. However, no check was 
made to see if the TPS was over design. The optimum panel configuration was deter- 
mined simply by varying panel size and looking either for minimum weight or cost of 
the TPS. The technique has been made more straight-forward in the new computer 
code by requiring that the program user determine the proper panel thicknesses in a 
series of computer runs by plotting, for example, the constraining temperature as a 
function of insulation thickness. In this manner, all the information obtained on TPS 
performance for the off-design thicknesses could be obtained, and in addition, the exact 
thickness to meet the temperature constraint can be determined. Observation of the 
stress values and design factors for panel and supporting substrate configurations 
allow reaching similar design conclusions for the stress supporting parts of the TPS 
configuration . 

For the case of the metallic panels, the thickness is still varied by increasing its 
previous value by a Newton-Raphson interpolation procedure to obtain all positive 
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values of the stress design factors. However, for the case of the RSI, although both 
discrete element and finite element stress analyses are performed, no automatic 
sizing is done - this is performed by the program user. 

Optimization and sensitivity studies are performed by the user by varying panel size, 
material properties, and configuration (six different metallic panel cross-section 
geometries are provided) in a series of computer runs. The program sizes panel 
and insulation thicknesses. An optimum design is then identified as the one giving 
either minimum weight or cost as a function of the parameters being varied for the 
investigation. Sensitivity studies are performed by noting the change in system weight 
or cost due to the variation in some independent variable such as trajectory or heating 
prediction method for an optimum panel configuration. 

Some typical results are available for an optimization study performed on the bottom 
centerline of the orbiter vehicle for an integrally stiffened panel composed of coated 
Columbium, Ren4 41, and Inconel 718. The panel is insulated with 6 pound per cubic 
foot Dynaflex. An underlaying aluminum tank constrained the system to a maximum 
temperature of 300° F here. The weight and manufacturing cost evaluated for a 
configuration that is 24 inches wide and of varying length is given in Figure 3-2. 

Cost of the coated Columbium panel comes out at approximately $700 per square foot 
and is relatively insensitive to panel length above a minimum value of approximately 
20 inches. For smaller panel lengths, the cost of the supporting structure becomes 
excessive and drives the total system cost up. The two additional materials (Inconel 
718 and Ren^41) indicated the same flat weight curve and relatively lower costs. 
Again, costs increase rapidly at the smaller panel sizes. The optimum value for 
both materials is a little over $200 per sq ft. 

3.4 SENSITIVITY STUDIES 

The prediction techniques used to design the TPS always have a degree of uncertainty 
caused by anything from the scatter of data used to correlate heat transfer prediction 
methods to the experimental error introduced while determining material properties. 
The system designer must account for these effects by either providing a margin of 
safety to allow for these contingencies or by limiting or refining mission constraints 
to avoid catastrophic failure due to the increased loads that may occur. The 
experience of previous TPS design work indicates that the compounded conservatism 
introduced by simply adding the margins required to provide system safety can lead 
to significant weight and performance penalties of the vehicle. To combat these 
penalties in both the weight and cost of the TPS, two tasks must be performed: (1) 
the penalties due to the performance prediction uncertainties must be established, 
and (2) a rational method of combining these penalties must be established since many 
of the uncertainties and associated penalties are not independent. It is the purpose of 
a sensitivity study to establish the penalties in system weight and cost due the 
uncertainties in system design. 
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Design uncertainties in the TPS itself can be separated into three categories: the 
environment, physical properties of the system itself, and operational and system 
characteristics. A partial list of these effects is given in Table 3-1. Parametric 
data gathered by establishing the weight and cost penalties due to variations in these 
uncertainties will establish the system sensitivities to these changes. 

Table 3-1. TPS Uncertainty Factors 


Environment 

Laminar heating 

Onset of transition to turbulent flow 

Onset of fully developed turbulent flow 

Turbulent heating 

Separated flow heating 

Local flow conditions 

Trajectory dispersions 

Roughness effects 

Vehicle aerodynamic attitude 

Local aerodynamic pressure loading 

Venting pressure 

Physical Properties 

Material properties, (including conductivity, density, specific heat, surface 
emissivity, etc) 

Material temperature limitations 
Mechanical properties 

Operations and System Characteristics 

Reentry trajectory corridor 

Attitude control system interactions 

Related systems and structure thermal limits 

Manufacturing, assembly and fabrication limitations 

Reuse constraints 

Inspection requirements 

Initial structural temperatures 
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SECTION 4 


CONCLUSIONS AND RECOMMENDATIONS 


1. The computer program developed under this contract provides the user with a 
powerful tool to determine the minimum weight/cost TPS for a given application 
from a variety of alternate TPS concepts and alternate materials for a given 
concept. 

2. An extensive parametric study should be performed on the current space 
shuttle vehicle to define the TPS concept, or combination of different concepts 
used at different locations on the vehicle which will yield the lowest cost TPS. 

3. Formal optimization procedures, utilizing automated iterative loops with their 
attendant curve fitting or extrapolation routines do not appear to be compatible 
with the requirements for a computer program with an acceptable cost. Both 
the cost of performing the program verification (assuming one is developed), 
and then the operating cost when applied to a problem appear to be excessive 
when compared with a man-in-the-loop procedure. 
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